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We study elastic systems such as interfaces or lattices, pinned by quenched disorder. To escape triviality 
as a result of "dimensional reduction", we use the functional renormalization group. Difficulties arise in the 
calculation of the renormalization group functions beyond 1-loop order. Even worse, observables such as the 
2-point correlation function exhibit the same problem already at 1-loop order. These difficulties are due to the 
non-analyticity of the renormalized disorder correlator at zero temperature, which is inherent to the physics 
beyond the Larkin length, characterized by many metastable states. As a result, 2-loop diagrams, which involve 
derivatives of the disorder correlator at the non-analytic point, are naively "ambiguous". We examine several 
routes out of this dilemma, which lead to a unique renormalizable field-theory at 2-loop order. It is also the 
only theory consistent with the potentiality of the problem. The /3-function differs from previous work and the 
one at depinning by novel "anomalous terms". For interfaces and random bond disorder we find a roughness 
exponent C = 0.20829804e + 0.006858e^ , e = 4 - d. For random field disorder we find C ^ e/3 and compute 
universal amplitudes to order O(e^). For periodic systems we evaluate the universal amplitude of the 2-point 
function. We also clarify the dependence of universal amplitudes on the boundary conditions at large scale. All 
predictions are in good agreement with numerical and exact results, and an improvement over one loop. Finally 
we calculate higher correlation functions, which turn out to be equivalent to those at depinning to leading order 



I. INTRODUCTION 



Elastic objects pinned by quenched disorder are central to the 
physics of disordered systems. In the last decades a consid- 
erable amount of research has been devoted to them. From 
the theory side they are among the simplest, but still quite 
non-trivial, models of glasses with complex energy landscape 
and many metastable states. They are related to a remarkably 
broad set of problems, from subsequences of randonipermu- 
tations in mathematics jll Ell, random matrices |5|l to 
growth models |6, 7, 's'.llJmTni [ll[ll[ll and Burgers 
turbulence in physics 1 15l llql . as well as directed polymers 
^ [TtIi and optimization problems such as sequence align- 
ment in biology fTHHIl^]. Foremost, they are very use- 
ful models for numerous experimental systems, each with its 
specific features in a variety of situations. Interfaces in mag- 
nets i2Tl I22I1 experience either short-range disorder (random 
bond RB), or long range (random field RF). Charge density 
waves (CDW) |23| or the Bragg glass in superconductors 
27, .28.1 are periodic objects pinned by disorder. 
The contact line of liquid helium meniscus on a rough sub- 
strate is governed by long range elasticity |29, 30, 31]. All 
these systems can be parameterized by a A^-component height 
or displacement field Ux, where x denotes the d-dimensional 
internal coordinate of the elastic object (we will use Uq to de- 
note Fourier components). An interface in the 3D random 
field Ising model has d = 2, = 1, a vortex lattice d = 3, 
N — 2, a contact-line d — 1 and = 1. The so-called di- 
rected polymer (d = 1) has been much studied J p2ll as it maps 
onto the Kardar-Parisi-Zhang growth model [6] for any N. 
The equilibrium problem is defined by the partition function 



Z = cxp(— 7i[u]/T) associated to the Hamiltonian 

nlu]^ J d'^x^iVu)^ + V{ux,x) , (1.1) 

which is the sum of an elastic energy which tends to suppress 
fluctuations away from the perfectly ordered state u = 0, 
and a random potential which enhances them. The resulting 
roughness exponent ( 



{{u{x) - u{x'))^) - Ix-xf ^ 



(1.2) 



is measured in experiments for systems at equilibrium (Ccq) 
or driven by a force /. Here and below (...) denote thermal 
averages and (...) disorder ones. In some cases, long-range 
elasticity appears e.g. for the contact line by integrating out 
the bulk-degrees of freedom |31], corresponding to \q\ 
in the elastic energy. As will become clear later, the random 
potential can without loss of generality be chosen Gaussian 
with second cumulant 



V{u, x)V{u', x') = R{u - u')5'^{x - x') 



(1.3) 



with various forms: Periodic systems are described by a pe- 
riodic function R{u), random bond disorder by a short-range 
function and random field disorder of variance a by R{u) ^ 
—cr\u\ at large u. Although this paper is devoted to equilib- 
rium statics, some comparison with dynamics will be made 
and it is thus useful to indicate the equation of motion 



rtdtUxt = cS/luxt + F{x, Uxt) + f 



(1.4) 



with friction ry. The pinning force is F{u, x) — ~duV{u, x) 
of correlator A{u) — —R"{u) in the bare model. 
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Despite some significant progress, the model has 
mostly resisted analytical treatment, and one often has to rely 
on numerics. Apart from the case of the directed polymer in 
1 + 1 dimensions (d = 1, N = JJ^here a set of exact and rig- 
orous results was obtained l2ll5L l33,. ,34. , 351. analytical meth- 
ods are scarce. Two main analytical methods exist at present, 
both interesting but also with severe limitations. The first one 
is the replica Gaussian Variational Method (GVM) |36]. It 
is a mean field method, which can be justified for N = oo 
and relies on spontaneous replica symmetry breaking (RSB) 
Although useful as an approximation, its validity 
at finite N remains unclear Indeed, it seems now generally 
accepted that RSB does not occur for low d and N. The re- 
maining so-called weak RSB in excitations may 
not be different from a more conventional droplet picture. An- 
other exactly solvable mean field limit is the directed polymer 
on the Cayley tree, which also mimics N ^ oo and there too 
it is not fully clear how to meaningfully expand around that 
limit |42, 43, 44|. The second main analytical method is the 
Functional Renormalization Group (FRG) which attempts a 
dimensional expansion around d = i llHlllEllillijl- The 
hope there is to include fluctuations, neglected in the mean 
field approaches. However, until now this method has only 
been developed to one loop, for good reasons, as we discuss 
below. Its consistency has never been checked or tested in any 
calculation beyond one loop (i.e. lowest order in e = 4 — d). 
Thus contrarily to pure interacting elastic systems (such as e.g. 
polymers) there is at present no quantitative method, such as 
a renormalizable field theory, which would allow to compute 
accurately all universal observables in these systems. 

The central reason for these difficulties is the existence 
of many metastable states (i.e. local extrema) in these sys- 
tems. Although qualitative arguments show that they arise 
beyond the Larkin length 14811 . these are hard to capture 
by conventional field theory methods. The best illustration 
of that is the so called dimensional reduction (DR) phe- 
nomenon, which renders naive perturbation theory useless 
mils EH El El in pinned elastic systems as well as 
in a wider class of disordered models (e.g. random field spin 
models). Indeed it is shown that to any order in the disorder 
at zero temperature T — Q, any physical observable is found 
to be identical to its (trivial) average in a Gaussian random 
force (Larkin) model, e.g. C — {A ~ d)/2 for RB disorder. 
Thus perturbation theory appears (naively) unable to help in 
situations where there are many metastable states. The two 
above mentioned methods (GVM and FRG) are presently the 
only known ways to escape dimensional reduction and to ob- 
tain non-trivial values for C (in two different limits but con- 
sistent when they can be compared |26t i28i |47J). The mean 
field method accounts for metastable states by RSB. This how- 
ever may go further than needed since it implies a large num- 
ber of pure states (i.e. low (free) energy states differing by 
0(T) in (free) energy). The other method, the FRG, captures 
metastability through a non-analytic action with a cusp sin- 
gularity. Both the RSB and the cusp arise dynamically, i.e. 
spontaneously, in the limits studied. 

The 1-loop FRG has had some success in describing pinned 
systems. It was noted by Fisher 146,1 within a Wilson scheme 



analysis of the interface problem in d = 4 — e that the coarse 
grained disorder correlator becomes non-analytic beyond the 
Larkin scale Lc, yielding large scale results distinct from 
naive perturbation theory. Within this approach an infinite set 
of operators becomes relevant in d < 4, parameterized by the 
second cumulant R{u) of the random potential. Explicit so- 
lution of the 1-loop FRG for R{u) gives several non-trivial 
attractive fixed points (FP) to 0(e) proposed in |461 to de- 
scribe RB, RF disorder and in |26, 28], periodic systems such 
as CDW or vortex lattices. All these fixed points exhibit a 
"cusp" singularity as R*"{u) — R*"{0) ~ \u\ at small 
The cusp was interpreted in terms of shocks in the renormal- 
ized force |54J, familiar from the study of Burgers turbulence 
(for d = 1, N = 1). The dynamical FRG was also devel- 
oped to one loop [55, 56, 57] to describe the depinning transi- 
tion. The mere existence of a non-zero critical threshold force 
fc ^ |A'(0+)| > is a direct consequence of the cusp (it 
vanishes for an analytic force correlator A(u)). Extension 
to non-zero temperature T > suggested that the cusp is 
rounded within a thermal boundary layer u ~ TL^^. This 
was interpreted to describe thermal activation and leads to a 
reasonable derivation of the celebrated creep law for activated 
motion ElEl- 

In standard critical phenomena a successful 1-loop calcu- 
lation usually quickly opens the way for higher loop compu- 
tations, allowing for accurate calculation of universal observ- 
ables and comparison with simulations and experiments, and 
eventually a proof of renormalizability. In the present context 
however, no such work has appeared in the last fifteen years 
since the initial proposal of 1.46,1 . a striking sign of the high 
difficulties which remain. Only recently a 2-loop calculation 
was performed ] 6^, '6I'] but since this study is confined to an 
analytic it only applies below the Larkin length and does 
not consistently address the true larg e scale critical behavior 
In fact doubts were even raised f47] about the validity of the 
e-expansion beyond order e. 

It is thus crucial to construct a renormalizable field the- 
ory, which describes statics and depinning of disordered elas- 
tic systems, and which allows for a systematic expansion in 
e = 4 — d. As long as this is not achieved, the physical mean- 
ing and validity of the 1-loop approximation does not stand 
on solid ground and thus, legitimately, may itself be called 
into question. Indeed, despite its successes, the 1-loop ap- 
proach has obvious weaknesses. One example is that the FRG 
flow equation for the equilibrium statics and for depinning are 
identical, while it is clear that these are two vastly different 
physical phenomena, depinning being irreversible. Also, the 
detailed mechanism by which the system escapes dimensional 
reduction in both cases is not really elucidated. Finally, there 
exists no convincing scheme to compute correlations, and in 
fact no calculation of higher than 2-point correlations has been 
performed. 

Another motivation to investigate the FRG is that it should 
apply to other disordered systems, such as random field 
spin models, where dimensional reduction also occurs and 
progress has been slow ]45, 62, 63, 64, 65]. Insight into model 
( II. H will thus certainly lead to progresses in a broader class 
of disordered systems. 



3 



In this paper we construct a renormalizable field theory for 
the statics of disordered elastic systems beyond one loop. The 
main difficulty is the non-analytic nature of the theory (i.e. 
of the fixed point effective action) at T = 0. This makes 
it a priori quite different from conventional field theories for 
pure systems. We find that the 2-loop diagrams are naively 
"ambiguous", i.e. it is not obvious how to assign a value to 
them. We want to emphasize that this difficulty already ex- 
ists at one loop, e.g. even the simplest one loop correction to 
the two point function is naively "ambiguous" . Thus it is not 
a mere curiosity but a fundamental problem with the theory, 
"swept under the rug" in all previous studies, but which be- 
comes unavoidable to confront at 2-loop order It originates 
from the metastability inherent in the problem. For the re- 
lated theory of the depinning transition, we have shown in 
companion papers |66, 67] how to surmount this problem and 
we constructed a 2-loop renormalizable field theory from first 
principles. There, all ambiguities are naturally lifted using 
the known exact property that the manifold only moves for- 
ward in the slowly moving steady state. Unfortunately in 
the statics there is no such helpful property and the ambi- 
guity problem is even more arduous. Here we examine the 
possible ways of curing these difficulties. We find that the 
natural physical requirements, i.e. that the theory should be 
(i) renormalizable (i.e. that a universal continuum limit exists 
independent of short-scale details), (ii) that the renormalized 
force should remain potential, and (iii) that no stronger singu- 
larity than the cusp in R"{u) should appear to two loop (i.e. 
no "supercusp "), are rather restrictive and constrain possible 
choices. We then propose a theory which satisfies all these 
physical requirements and is consistent to two loops. The 
resulting /3-function differs from the one derived in previous 
studies 1 60, 61] by novel static "anomalous terms". These are 
different from the dynamical "anomalous terms" obtained in 
166^ .67, .681 showing that indeed depinning and statics differ 
at two loop, fulfilling another physical requirement. 

We then study the fixed points describing several universal- 
ity classes, i.e. the interface with RB and RF disorder, the ran- 
dom periodic problem, and the case of LR elasticity. We ob- 
tain the O(e^) corrections to several universal quantities. The 
prediction for the roughness exponent for random bond dis- 
order has the correct sign and order of magnitude to notably 
improve the precision as compared to numerics in d — 3, 2 
and to match the exact result C, ~ 2/3 in d = 1. For ran- 
dom field disorder we find = e/3 which, for equihbrium 
is likely to hold to all orders. By contrast, non-trivial correc- 
tions of order O(e^) were found for depinning |66, 6*^. The 
amplitude, which in that case is a universal function of the 
random field strength is computed and it is found that the 2- 
loop result also improves the agreement as compared to the 
exact result known |69] for d = {). For the periodic CDW 
case we compare with the numerical simulations in d = 3 and 
obtain reasonable agreement. Some of the results of this paper 
were briefly described in a short version |66] and agree with a 
companion study using exact RG [70., 71]). 

Since the physical results also seem to favor this theory we 
then look for better methods to justify the various assump- 
tions. We found several methods which allow to lift ambigu- 



ities and all yield consistent answers. A detailed discussion 
of these methods is given. In particular we find that correla- 
tion functions can be unambiguously defined in the limit of 
a small background field which splits apart quasi-degenerate 
states when they occur. This is very similar to what was found 
in a related study where we obtained the exact solution of the 
FRG in the large N limit f72]. Finally, the methods introduced 
here will be used and developed further to obtain a renormaliz- 
able theory to three loops, and compute its /?-function in JtsIi . 
Let us mention that a first principles method which avoids am- 
biguities is to study the system at T > 0. However, this turns 
out to be highly involved. It is attempted via exact RG in 
and studied more recently in |74, 75] where a field theory of 
thermal droplet excitation was constructed. A short account 
of our work has appeared in Il66ll . and a short pedagogical in- 
troduction is given in L76.1 . 

The outline of this paper is as follows. In Section IIIl we 
explain in a detailed and pedagogical way the perturbation 
theory and the power counting. In Section |III] we compute 
the 1-loop (Section IlII A> and 2-loop (Section llIIB> correc- 
tions to the disorder The calculation of the repeated 1-loop 
counter-term is given in Section IIII CI In Section IIII Dl we 
identify the values for ambiguous graphs. This yields a renor- 
malizable theory with a finite /3-function, which is potential 
and free of a supercusp. The more systematic discussion of 
these ambiguities is postponed to Section|V] We derive the j3- 
function and in Section HVl present physical results, exponents 
and universal amplitudes to O(e^). Some of these quantities 
are new, and have not yet been tested numerically. In Section 
|V]we enumerate all the methods which aim at lifting ambigu- 
ities and explain in details several of them which gave con- 
sistent results. In Section fVIl we detail the proper definition 
and calculation of correlation functions. In Appendix IaI and 
|n]we present two methods which seem promising but do not 
work, in order to illustrate the difficulties of the problem. In 
Appendix|F|we present a summary of all one and 2-loop cor- 
rections including finite temperature. In AppendixiDlwe give 
details of calculations for what we call the sloop elimination 
method. 

The reader interested in the results can skip Section IIII and 
Section|III]and go directly to Section HVI The reader interested 
in the detailed discussion of the problems arising in this field 
theory should read Section [VI 



II. MODEL AND PERTURBATION THEORY 

A. Replicated action and effective action 

We study the static equilibrium problem using replicas, i.e. 
consider the partition sum in presence of sources: 



m = j n^^K] exp \^S[u] + j E-?>"j ' (2.1) 



from which all static observables can be obtained. The action 
S and replicated Hamiltonian corresponding to are 

n[u] _ _i_ 

^ 2T 



S[u] 



T 



a\2 



[(v<) 



ab 



(2.2) 



a runs from 1 to n and the limit of zero number of replicas n = 
is implicit everywhere. We have added a small mass which 
confines the interface inside a quadratic well, and provides 
an infrared cutoff. We are interested in the large scale limit 
m — > 0. We will denote 



= / d^a 



(2.3) 



(2.4) 



For periodic systems the integration is over the first Brillouin 
zone. A short-scale UV cutoff is implied at (7 ^ A, but for 
actual calculations we find it more convenient to use dimen- 
sional regularization. We also consider the effective action 
functional T[u\ associated to S. It is, as we recall f77ll73l . the 
Legendre transform of the generating function of connected 
correlations W[j] = In thus defined by eliminating j in 
T[u]=ju-W[i\,W'[j]^u. 

If we had chosen non-Gaussian disorder additional terms 
with free sums over p replicas (called p-replica terms) corre- 
sponding to higher cumulants of disorder would be present in 
(12. 2> . together with a factor of l/T^. These terms are gen- 
erated in the perturbation expansion, i.e. they are present in 
T[u]. We do not include them in M.2\ because, as we will see 
below, these higher disorder cumulants are not relevant within 
(conventional) power counting, so for now we ignore them. 
The temperature T appears explicitly in the replicated action 
(12. 2> . although we will focus on the T = limit. 

Because the disorder distribution is translation invariant, 
the disorder term in the above action is invariant under the 
so called statistical tilt symmetry OH (STS), i.e. the shift 
u% u% + (jx- One implication of STS is that the 1-replica 
replica part of the action (i.e. the first line of is uncor- 
rected by disorder, i.e. it is the same in T[u\ and S[u\ isoll . 
Since the elastic coefficient is not renormalized, we have set 
it to unity. 

B. Diagrammatics, definitions 

We first study perturbation theory, its graphical representation 
and power counting. Everywhere in the paper we denote the 
exact 2-point correlation by Cab{x — y), i.e. in Fourier: 

{u-gul,) = {27:YS''iq + q')Cabiq) (2.5) 

while the free correlation function (from the elastic term) used 
for perturbation theory in the disorder is denoted by Gab{x — 
y) — 5abG{x — y) and reads in Fourier: 

(2.6) 
(2.7) 



{uy^,)^^{2^Y5%q + <^)Gab{q) 
T 



Gab{q) 



q^ + m? 



Sab 



FIG. 1: Each diagram with unsplitted vertices contains several dia- 
grams with splitted vertices: here the 1-loop unsplitted diagram (top) 
generates three possible topologically distinct splitted diagrams, two 
(shown here, bottom) are 2-replica terms, the third one, i.e. (a) in Fig. 
J2j is a three replica term 



which is represented graphically by a line: 

T6ab 



(2.8) 

q- -t- riL- 

Each propagator thus carries one factor of G{q) = T / {q^ + 
m?). Each disorder interaction vertex comes with a factor of 
l/T^ and gives one momentum conservation rule. Since each 
disorder vertex is a function, an arbitrary number of lines can 
come out of it. k lines coming out of a vertex result in k 
derivatives R^'^^ after Wick contractions 



(2.9) 



Since each disorder vertex contains two replicas it is some- 
times convenient to use "splitted vertices" rather than "un- 
splitted ones". Thus we call "vertex" an unsplitted vertex and 
we call a "point" the half of a vertex. 



E 



R{Ua - Ub) 

2T2 



(2.10) 



Each unsplitted diagram thus gives rise to several splitted dia- 
grams, as illustrated in Fig.^ 

One can define the number of connected components in a 
graph with splitted vertices. Since each propagator identifies 
two replicas, a p-replica term contains p connected compo- 
nents. When the 2-points of a vertex are connected, this ver- 
tex is said to be "saturated". It gives a derivative evaluated at 
zero i?('^')(0). Standard momentum loops are loops with re- 
spect to unsplitted vertices, while we call "sloops" the loops 
with respect to points (in splitted diagrams). This is illustrated 
in Fig. (|2} The momentum 1-loop and 2-loop diagrams which 
correct the disorder at T = are shown in Fig. |3] (unsplitted 
vertices). There are three types of 2-loop graphs A, B and C. 
Since they have two vertices (a factor R/T^ each) and three 
propagators (a factor of T each) the graphs E and F lead to 
corrections to R proportional to temperature and will not be 
studied here (see however Appendix|F|i. 

It is important to distinguish between fully saturated dia- 
grams and functional diagrams. The FS diagrams are those 
needed for a full average, e.g. a correlation function. There all 
fields are contracted and one is only left with the space depen- 
dence. These are the standard diagrams in more conventional 
polynomial field theories such as 0"*. Then all vertices are 
evaluated at it = 0, yielding products of derivatives i?(*'''(0). 
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(a) 



(c) 



(b) 




tain saturated vertices, whose space and field dependence dis- 
appears (such as (c) in Fig.O and that the Hmit n ^ does 
not produce constraints. An example is the calculation of r[M] 
since one can always attach additional external legs to any 
point by taking a derivative with respect to the field u. 



FIG. 2: Graphs (a) (a 1-loop diagram) and (b) (a 2-loop diagram) 
each contains three connected components. Since each contain one 
"sloop" they are both three replica terms proportional to T. The left 
vertex on diagram (c) is "saturated" : replica indices are constrained 
to be equal and thus the diagram does not depend on the left space 
point. 




C 



OO CO 
<2> O - 



FIG. 3: unsplitted diagrams to one loop D, one loop with inserted 
1-loop counter- term G and 2-loop diagrams A, B, C, E and F. 



These are also the graphs which come in the standard expan- 
sion of r[u] in powers of u which generate the "proper" or 
"renormalized" vertices, i.e. the sum over all 1 -particle irre- 
ducible graphs with some external legs, from which all corre- 
lations can be obtained. Note that in the fully saturated dia- 
grams there can be no free point, all points in a vertex have 
to be connected to some propagator (and to some external 
replica) otherwise there is a free replica sum yielding a fac- 
tor of n and a vanishing contribution in the limit of n = 0. 

However, since we have to deal with a function B.{u) we 
will more often consider functional diagrams. A functional 
diagram still depends on the field u. It can depend on u at 
several points in space (multi-local term), as for example: 



E 

abc 



abc 

(2.11) 

Such a graph with p connected components corresponds to a 
p replica functional term. Or it can represent the projection of 
such a term onto a local part, as arises in the standard operator 
product expansion (OPE): 



abc 



ui) R'{u% - u%) 



2T2 



2T2 



T / G{x^y) 
(2.12) 

Typically using functional diagrams we want to compute the 
effective action functional T[u], or its local part, i.e. its value 
for a spatially uniform mode u° = u", which includes the 
corrections to disorder Specifying the two replicas on each 
connected component, one example of a 1 -particle irreducible 
diagram producing corrections to disorder is 



rp4 



R" {u'' ^ u'')R" [u" 



Jq 

(2.13) 

The complete analysis of these corrections will be made in 
Section JIII> . Finally, note that functional diagrams may con- 



C. Dimensional reduction 

If we consider fully saturated diagrams and analytic R{u) we 
find trivial results. This is because at T = the model exhibits 
the property of dimensional reduction i2lll49ll5(il5lll5^l53ll 
(DR) both in the statics and dynamics. Its "naive" pertur- 
bation theory, obtained by taking for the disorder correlator 
R{u) an analytic function of u has a triviality property. As is 
easy to show using the above diagrammatic rules (see a typ- 
ical cancellation due to the "mounting" construction in Fig. 
13 see also Appendix D in Ref. |70]) the perturbative expan- 
sion of any correlation function (W^u°^^,)s (of any analytic 
observable) in the derivatives i?*^*^) (0) yields to all orders the 
same result as that obtained from the Gaussian theory set- 
ting R{u) = R"{Q)v? /2 (the so called Larkin random force 
model). The 2-point function thus reads to all orders: 



-i?"(0) 



(2.14) 



(at T = correlations are independent of the repUca indices 
tti). This dimensional reduction results in a roughness expo- 
nent C = (4 — c?)/2 which is well known to be incorrect. One 
physical reason is that this T = perturbation theory amounts 
to solving in perturbation the zero force equation 



{-V'^ +m'^)u + F[x,u) = 



(2.15) 



This, whenever more than one solution exists (which cer- 
tainly happens for small m) is cl early not identical to finding 
the lowest energy configuration lll02ll . Curing this problem 
within the field theory, is highly non-trivial. Coarse graining 
within the FRG up to a scale at which the renormalized dis- 
order correlator R{u) becomes non-analytic (which includes 
some of the physics of multiple extrema) is one possible route, 
although understanding exactly how this cures the problem 
within the field theory is a difficult open problem. 

It is important to note that dimensional reduction is not the 
end of perturbation theory, since saturated diagrams remain 
non-trivial at finite temperature, so one way out is to study 
T > 0. This is not the route chosen here, instead we will 
attempt to work at T = with a non-analytic action and focus 
on functional diagrams which remain non-trivial. 

D. Power counting 

Let us now consider power counting. Let us recall the con- 
ventional analysis within e.g. the Wilson scheme |46, 47^. 
The elastic term is invariant under x bx, u h'^u and 
T b^T, with 6 = d-2 + 2C. (is for now undetermined. 
Under this transformation the disorder function R is multi- 
plied by fo'^^^e _ ^4-d+2C jj becomes relevant for d < 4, 
provided ( < {4 — d)/2 which is physically expected (for 
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(a) (b) (c) 




a 9 « 

• • h 

(d) 

FIG. 4: Calculation of the 2-point function for analytic R(u). Due to 
DR only the first diagram (a) survives. Diagrams (b) and (c) cancel 
because by shifting the line one gets a minus sign. Diagram (c) is 
proportional to i?"'(0)^ and vanishes in an analytic theory. Similar 
cancellations occur to all orders. 

instance in the random periodic case, = is the only pos- 
sible choice, and for other cases ^ — 0(e)). The rescaled 
dimensionless temperature term scales as —mdmT — —9T 
(see below) and is formally iiTelevant near four dimension. In 
the end C, will be fixed by the disorder distribution at the fixed 
point. 

To be more precise, we want to determine in the field theo- 
retic framework the necessary counter-terms to render the the- 
ory UV finite as d ^ 4. The study of superficial divergences 
usually involves examining the irreducible vertex functions 
(IVF): 

r„...„fe) = n^rM|_^ (2.16) 

with Eu external fields u (at momenta qi, i = 1, --Ey). The 
perturbation expansion of a given IVF to any given order in 
the disorder is represented by a set of 1 -particle irreducible 
(IPI) graphs (in unsplitted diagrammatics). Being the deriva- 
tive of the effective action they are the important physical ob- 
jects since all averages of products of fields u can be expressed 
as tree diagrams of the IVF. Finiteness of the IVF thus implies 
finiteness of all such averages. 

However since V[u] is non-analytic in some directions (e.g. 
for a uniform mode u% = m°), derivatives such as ( 12.161 1 may 
not exist at g = 0, and we have to be more general and con- 
sider /Mncf/onaZ diagrams. The (disorder part of the) effective 
action is the sum of fc-replica terms, noted [u] 

T[u] = Y,T-''Tk[u]. (2.17) 

k>2 

Each Ffc [u] is the sum over IPI graphs with k connected com- 
ponents (using splitted vertices), and itself depends on T as 

rfeM^^T'rfe^M, (2.18) 

where I is the number of sloops. Thus at T = there are 
no sloops and Tk[u] = r/j.;=o[u] is the sum over IPI tree 
graphs with k connected components (trees in replica-space, 
not position-space). 



Let us compute the superficial degree of UV divergence 5 
of a functional graph entering the expansion of the local part 
of the effective action. We denote v the number of unsplitted 
disorder vertices, / the number of internal lines (propagators), 
L the number of loops and I the number of sloops. One has 
the relations 

2v + l = k + I (2.19) 
v + L = l + I . (2.20) 

The total factors of T are T^-^" = T'-''. AtT = 0(1 ^ 0) 
the superficial degree of UV divergence is thus 

5 = dL-2I = d- k{d-2) + {d- 4)v . (2.21) 

Thus in d = 4 the only graphs with positive superficial degree 
of divergence are for k — 1 (quadratic ^ A^), and k — 2 (log 
divergence), k = 1 corresponds to a constant in the free en- 
ergy. Because of STS all single replica terms are uncorrected 
and there is no wave-function renormalization in this model. 

Thus to renormalize the T = theory we need a priori 
to look only at graphs with p = 2 connected components, 
which by definition are those coiTecting the second cumulant 
R{u), compute their divergent parts, and construct the proper 
counter-term to the function R{u). As mentioned above, 
higher cumulants are irrelevant by power counting, and are 
superficially UV-finite. The graphs which contribute to the 2- 
replica part T2[u] have L loops with L = 1 + v + I. At zero 
temperature, I = 0, thus L = 1 + v. The loop expansion thus 
coiTesponds to the expansion in power of R{u) and, as we will 
see below, to an e-expansion. More generally using the above 
relation one has, schematically 

L>max(l,2+;-fe) 

(2.22) 

where the number of internal lines gives the total number of 
derivatives acting on an argument u of the functions R. For 
instance, the 2-replica part at T = is a sum over L-loop 
graphs of the type 

L>1 

If one now considers T > one finds that 6 — d — k{d ~ 
2) + (rf - 4)u + (d - 2)1. Each additional power of T yields 
an additional quadratic divergence, more generally a factor of 
2i^(i-2 -pjjyg j-Q obtain a theory where observables are finite 
as A ^ oo one must start from a model where the initial 
temperature scales with the UV cutoff as 

T = fk^-'^ . (2.24) 

This is similar to 0'*-theory where it is known that a (f)^ term 
can be present and yields a finite UV limit (i.e. does not spoil 
renormalizability) only if it has the form ggcf)^ / A"^^^ . Such a 
term, with precisely this cutoff dependence, is in fact usually 
present in the starting bare model, e.g. in lattice spin models. 
It then produces only a finite shift to 54 without changing uni- 
versal properties Iil03i] . Here each factor of T comes with a 
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(a) 



(b) 



FIG. 5: The two 1-loop diagrams with splitted vertices and the cor- 
responding diagram with standard (i.e. unsplitted) vertices. 



factor of A^"'' which compensates the UV divergence from 
the graph. Thus the finite-T theory may also be renormaHz- 
able. Computing the resulting shift in R{u) to order by 
resumming the diagrams E and F of Fig. |3] and all similar 
diagrams to any number of loops has not been attempted here 
(see however AppendixIQ. The "finite shift" here is, however, 
much less innocuous than in (^"'-theory since it smoothes the 
cusp . The effects of a non-zero temperature are explored in 

ilTilzilzlil. 

One can use the freedom to rescale u by m^^. The dimen- 
sionless temperature T = Tm^ is then defined. The disor- 
der term in r[u] is then is as in i2.2i with R{u) replaced by 
rrf~'^^ R{u'nrfi) in terms of a dimensionless rescaled function 
J? of a dimensionless rescaled argument. This will be further 
discussed below. 



III. RENORMALIZATION PROGRAM 

In this section we compute the effective action to 2-loop order 
at T = 0. We are only interested in the part which contains 
UV divergences as d ^ 4. We know from the analysis of 
the last section that we only need to consider the local k — 2 
2-replica part, i.e. the corrections to R{u). These L — 1 and 
L = 2 loop corrections contain v ^ L + 1 vertices. Higher v 
yields higher number of replicas. 

A. 1-loop corrections to disorder 

To one loop at T = there is only one unsplitted diagram 
V = 2, corresponding to two splitted diagrams (a) and (b) as 
indicated in figure|5] Both come with a combinatorial factor 
of 1/2! from Taylor-expanding the exponential function and 
1 /2 from the action, (a) has a combinatoric factor of 2 and (b) 
of 4. Together, they add up to the 1-loop correction to disorder 



S^R{u) 

ir = h 



^R"iuf - R"iO)R"{u) 
1 



(3.1) 



g (q^ + m^)^ 
1 



r 2 



tT 2 



= O 



(3.2) 



Note that (b) has a saturated vertex, hence the factor i?"(0). 
This does not lead to ambiguities in the 1-loop /3-function, 
since the FRG to one loop yields a discontinuity only in the 
third derivative and R"{u) remains continuous. 



B. 2-loop corrections to disorder 

There are only three graphs correcting disorder at T = with 
L = 2 loops and w = 3 vertices. They are denoted A, B and 
C and we will examine each of them. 
We begin our analysis with class A. 

1 . Class A 

The possible diagrams with splitted vertices of type A are di- 
agrams (a) to (f) given in Fig. The resulting correction to 
R{u) is written as: 



8'^R{u) 



1 2 
3! 23 



3(23)^(a + 

b + c + d 



f) 



(3.3) 



where the combinatorial factors are: 1/3! from the Taylor- 
expansion of the exponential function, 2/2^ from the explicit 
factors of 1/2 in the interaction, a factor of 3 to chose the 
vertex at the top of the hat, and a factor of 2 for the possible 
two choices in each of the vertices. Furthermore below some 
additional combinatorial factors are given: A factor of 2 for 
generic graphs and 1 if it has the mirror symmetry with respect 
to the vertical axis. Each diagram symbol denotes the diagram 
including the symmetry factor The first two graphs are: 



a = -R"{{))R"'{uflA 
b = R"{u)R"'{uflA ■ 



(3.4) 
(3.5) 



To obtain the sign one can choose an "orientation" in each ver- 
tex {ua — Ub), the final result does not depend on the choice. 
The minus sign in a comes because the two legs enter on op- 
posite points in the top vertex. Define the 2-loop momentum 
integral (see Appendix A in Ref. j^]) 



Ia ■■= 



1 



1 



1 



qi ■Jq2 



2^2 



(3.6) 



Graphs a and h are non-ambiguous. They are the only contri- 
butions in an analytic theory. The other graphs are 

c = 2\^R"'{Q)R"{Q)R!"{u)Ia (3.7) 

d = 2XdR'"{0)R"{u)R"'{u)lA (3.8) 

e=-Ae(i?"'(0+))2i?"(u)/^ (3.9) 

f = 2XfR"'{0fR"{u)lA (3.10) 

and vanish if R{u) is analytic (since then i?"'(0) = 0) but a 
priori should be considered when R{u) is non-analytic. We 




A 



B 




FIG. 6: The 3 possible 2-loop unsplitted graphs correcting disorder 

atr = o 
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FIG. 7: Graphs at 2-loop order in the form of a hat (class A in figure 
IIIIB 1> contributing to 2-replica terms. 
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h 




J 



FIG. 8: 2-Ioop diagrams of class B 



have indicated their "natural" sign and amplitude (e.g. sym- 
metry factor setting = 1) but have introduced factors to 
recall that they are ambiguous: since i?"'(0+) = —R"'{0^) 
one is confronted to a choice each time one saturates a vertex 
and there is no obvious way to choose the sign at this stage. 
We recall that we have defined saturated vertices as vertices 
evaluated at u = while unsaturated vertices still contain u 
and do not lead to ambiguities. 

At this stage we will not discuss in detail how to give a 
definite values to these contributions to disorder. This will 
be done in Section IVl We will just use the most reasonable 
assumptions, which will be reevaluated, and justified later. A 
natural step is to set 



c = d = 0, 



(3.11) 



since these graphs cannot coiTect R{u) as they are odd func- 
tions of u, which yields no contribution when inserted into the 
action Y.ab ^("^^ ~ ut). 



2. Class B 

We now turn to graphs of type B (bubble-diagrams), g to I 
represented in Fig.|8] We use the same convention as in ( I3.3l l. 
and start with the combinatorics. There are 3 ways to choose 
the vertex in the middle. Upon splitting the vertices, for i and 
j there are only two choices at the middle vertex whereas for g 
there are four choices. There are also four choices for h, k and 
I. There, one must also choose the rightmost vertex, leading 






FIG. 9: 2-loop diagrams of class C 
to an extra factor of 2. The final result is 

h = -R"{u)R""{u)R"iO)lf 

^^J^\R"'\u)R"{OfI! 

k = -XkR"iu)R"{0)R""iO)lf 
I = XiR"{u)R"{0)R""{0)lf . 



(3.12) 

(3.13) 

(3.14) 

(3.15) 
(3.16) 



Only k and I are ambiguous but it is also natural to set: 

k + l = 0, (3.17) 
which we do for now, and discuss later. 

3. Class C 

Diagrams m, n, p, q of class C are represented in Fig.|9] 



m = ci\nR"iO)R""{0)R"{u)ItlT 
n = -ciA„i?"(0)i?""(0)i?"(M)/t/T 
p = c2\R!"'{Q)R!'{ufltlT 



q 



with 



= -C2XgR""iO)R"{u)^ItlT 



It = 



It 



1 



There it is natural to assume 



m + n = 
p + q = 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



(3.22) 
(3.23) 



(3.24) 
(3.25) 



which we do for now and discuss it later. This leaves no cor- 
rection to disorder from graphs C, as is the case for depinning 
i67il . This is fortunate, since the integral It has a quadratic 
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UV-divergence in d = 4, while It is UV-finite. Physically, it 
is unlikely that these could enter physical observables as the 
tadpole divergence can usually be eliminated by proper field 
reordering (normal-ordering) or vacuum subtraction. 

To summarize, for the equilibrium statics at T = in per- 
turbation of i? = R{u), the contributions to the disorder to 
one and two loops, i.e. the corresponding terms in the effec- 
tive action r[u, ti] are 



\R"{uf-R"{m"{u) 



h 



(3.26) 



(5'i?(u) = [R!"{uf{R!'[u) - i?"(0))] Ia 

+ \[{R"{u)-R"{Q)YR""{u)]iI 
-\R"\Q+fR"{u)lA . (3.27) 

We have allowed for a yet undetermined constant A = Ae — 
2A/. We now show that requiring renormalizability allows to 
fix A. 

C. Renormalization method to two loops and 
calculation of counter-terms 

Let us now recall the method, also used in our study of depin- 
ning 1 67], to renormalize a theory where the interaction is not 
a single coupling-constant, but a whole function, the disorder- 
correlator R{u). We denote by i?o the bare disorder - this is 
the object in which perturbation theory is carried out - i.e. one 
considers the bare action M.2\ with R Rq. We denote here 
by R the renormalized dimensionless disorder i.e. the corre- 
sponding term in the effective action T[u] is ni'^R (i.e. the 
local 2-replica part of r[u]). Symbolically, we can write 



S[u] Rq 
T[u\ ^ m'^R . 



(3.28) 
(3.29) 



We define the dimensionless symmetric bilinear 1-loop and 
trilinear 2-loop functions (see ( I3.26> and ( I3.27t ) such that 



s'^^\r,r) 
s^^\r,r,r) 



m'S^R 
m'6^R 



(3.30) 
(3.31) 



They can be extended to non-equal argument using f{x,y) :— 
\ [fi.x + y,x + y)- f{x, x) - f{y, y)] and a similar expres- 
sion for the trilinear function. Whenever possible we will use 
the shorthand notation 5^^^\R) = S'-^^R,R) and S'-^^R) ^ 
R). The expression of R obtained perturbatively in 
powers of Rq at 2-loop order reads: 

R = 77i-'i?o + 6'-^\m-'RQ) + 6'-^\m-'RQ) + 0{R^) . 

(3.32) 

It contains terms of order 1/e and This is sufficient to 

calculate the RG-functions at this order In principle, one has 
to keep the finite part of the 1-loop terms, but we will work in 
a scheme, where these terms are exactly 0, by normalizing all 
diagrams by the 1-loop diagram. Inverting ( I3.32> yields 



R - S'-^^ (R) - (5(2) (i?) + ^(i^i) (i?) + 



where S^'^'^^R) is the 1-loop repeated counter-term: 

(i?) = 2^(1' {R, S'^^l [R, R)) . (3.34) 

The /3-function is by definition the derivative of R at fixed Rq. 
It reads 



md„iR\ =e 



m-''Ro + 2S^^Hm-'Ro) 
+3(5(2)(m-<^i?o) + -.-j ■ (3.35) 



Using the inversion formula (I3.33> . the /3-function can be 
written in terms of the renormalized disorder R: 



R + 6^^\R) 
+2,5(2)(i?)-5(i^i)(i?) 



(3.36) 



In order to proceed, let us calculate the repeated 1-loop 
counter-term 6^'^{R). We start from the 1-loop counter-term 
( I3.26> . which has the bilinear form 

-r{u)g"{0)]h (3.37) 
with the dimensionless integral L := lA , ; we will use the 

^ \m—l 

same convention for Ia := Thus S^'^{R) reads 



= -i?"(0))i?"'(M)2 

+ iR"{u) - R"{{))fR!"'{u) - R"\0+fR"{u) 



It- 
(3.38) 



In the course of the calculation the only possible ambiguity 
could come from 



1 



5"(0)= ^-R"{uf-R"{0)R"{u) 



R"\uf - R""{u){R"{u) - i?"(0)) 
R"'{Q+f 



(3.39) 



(3.33) 



but there is no ambiguity since the function R"'[uf' is con- 
tinuous at u with value R"'{Q+f = R"'{0- f. This is 
exactly the same calculation as is done to one loop when com- 
puting the non-trivial fixed point for the pinning force correla- 
tor A(u) = -R"{u) yielding = (e - 2C)A(0) - A'(0+)2. 
Thus there is no doubt that the graph G with the 1-loop 
counter-term inserted in a 1-loop diagram is non-ambiguous. 



10 



D. Final /3-function, renormalizability and 
potentiality 

The 2-loop /3-function ( I3.36> then becomes with the help of 



mdmR{u) — eR{u) 



(eii) 



-R"iuf - R"{0)R"{u) 



- [iR"{u) - R"{0))R"'{uf] e{2h - if 
'{R"'{0+)fR"iu)e(2XiA- if) ■ (3.40) 



The first result is that, apart from the last "anomalous" term, 
the 1/e^-terms cancel in the corrections to disorder In the 
terms coming from graphs A this works because, as we re- 
call, Ia = + 472- + 0(e^))(e/i)2 so that the combination 
e{2lA — if ) is finite. Graphs B cancel completely since we 
have chosen as counter-term the full 1-loop graph. So for an 
analytic theory the above /^-function would be finite. This 
however is incomplete, since the flow of such a /3-function 
leads to a non-analytic R{u) above the Larkin scale. 

Thus we must consider the last, "anomalous" term in ( I3.40> . 
It clearly appears that the only value of A compatible with the 
cancellation of the poles is 



A ^ 1 



(3.41) 



leading to a finite /3-function. Thus the requirement that the 
theory be renormalizable (i.e. yield universal large scale re- 
sults independent of the short-scale details) fixes the value 
A = 1. Note that the cancellation of the graphs B also works 
thanks to ( 13.17k 

It is interesting to compare with what happens at depinning. 
There the cancellation of the 1 / e^-terms in the anomalous part 
is more complicated but automatic. It requires a consistent 
evaluation of all anomalous non-analytic diagrams. In the 
depinning theory the cancellation was unusual: a non-trivial 
bubble diagram (called 13 in ll67ll ') was crucial in achieving the 
cancellation. In the statics the 2-loop bubble diagrams of type 
B appear to be simply the square of the 1-loop ones which 
is the usual situation. This however is clearly a consequence 
of (|^^} so the previous experience with depinning indicates 
that care is required and we will discuss some justification for 
(I3.17> below. 

In the search for a fixed point it is convenient to write the 
/3-function for the rescaled function R{u) defined through 



R{u) = -^m-'^'^ R{um'^) 
eh 



(3.42) 



which amounts to rescale the fields u by m'^. Note that this 
is a simple field rescaling and different from standard wave- 
function renormalization, since as mentioned above there is 
none in this theory due to STS. We have also included the 
1-loop integral factor to simplify notations and further calcu- 
lations (equivalently it can be absorbed in the normalization 
of momentum or space integrals). With this, the /3-function 



takes the simple form: 

- mdmR{u) = (e - 4C)^(m) + (uR'{u) 

^R"{uf - i?"(0)i?"(u) 



{R"{u) - R"{Oj)R"'{uf 



A = 1 



-^X{R'"{0+jfR"{u) 
X 1 . 



(3.43) 
(3.44) 



We have left a A for future use, but its value in the theory we 
study here is set to 1. Also for convenience we have intro- 
duced 



X = 



2e{2lA~lf) 
{eh? 



(3.45) 



which is X = 1 + 0(e) in the e-expansion studied here, but 
has a different value for LR elasticity, see below. In fact it 
is shown in Appendix lEI that lim^^o X is independent of the 
particular infrared cutoff procedure (here a massive scheme). 
Although the global rescaling factor of R, e/i, has 0(e) cor- 
rections which depend on the infrared cutoff chosen, the FRG 
equation above does not depend on it. Note that the above 
equation remains true in fixed dimension, with the appropri- 
ate value for X, up to terms of order R^. 

We will see that the value A = 1 in ( I3.43> has other highly 
desirable properties. First this value is the only one which 
guarantees that the non-analyticity in R{u) does not become 
more severe at two loops than it is at one loop. Let us take one 
derivative of ( I3.43> and take u 0+. One finds: 

-ma„,i?'(0+) = (e - 3C)i?'(0+) + ^(1 - \)R"' {Q+f . 

(3.46) 

Thus if A 7^ 1 the cusp in R" and the resulting finite value of 
R'" (0+) immediately creates a cusp in R! . The singularity has 
become worse! We call this a supercusp. It must be avoided 
in the statics (see also discussion in SectionlVl. Interestingly 
it does occur in the driven dynamics, where it is a physical 
signature of irreversibility. 

Indeed this property is intimately related to another highly 
desirable property of the statics: potentiality. This prop- 
erty is more conveniently described by considering the flow 
equation for the (rescaled) correlator of the pinning force 
A(u) — —R"{u), the second derivative of ( 13.431 1: 

- m9,„A(u) = (e - 2C)A(u) + CuA'{u) 
1 



2 
1 

'2 

4(A'(0+))2A"(u) 



(A(u)-A(0))A'(u)2 



(3.47) 



Formally, this equation could have been obtained directly 
from a study of the dynamical field theory. Such an equation 
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was indeed obtained at depinning but with a different value of 
A: 



A, 



dcp 



-1 



(3.48) 



which shows that statics and dynamics differ not at one, but 
at two loops. Integrating the equation for A(u) once yields 
a non-zero fixed point value for j A(u) unless A = 1. Po- 
tentiality on the other hand requires that the force remains the 
derivative of a potential and that, for short-range disorder (e.g. 
RB for interface) one must have J A(u) = 0. While violating 
potentiality is desirable at depinning where irreversibility is 
expected, this would be physically incorrect in the statics, and 
thus again points to the value A = 1 as the physically correct 
one. 

Thus we will for now assume that this is the correct theory 
of the statics and explore its consequences in the next section. 
In section|y]we will provide better justifications, and explain 
our understanding of the tantalizing problem of ambiguous di- 
agrammatics in the non-analytic theory of pinned disordered 
systems. Especially we will present methods, which satisfy 
all the above constraints of renormalizability, absence of a su- 
percusp and potentiality up to 3-loop order I 73ll . 

IV. ANALYSIS OF FIXED POINTS AND 
PHYSICAL RESULTS 

The FRG-equation derived above describes several different 
physical situations, and admits a small number of fixed-point 
functions R*{u) describing a few universality classes. The 
fixed point associated to a periodic disorder correlator de- 
scribes single component periodic systems (such as charge 
density waves). The fixed point associated to a short-range 
(exponentially decaying) correlator R{u) describes a class of 
systems with so called random bond disorder. There is also 
a family of fixed points associated to long range, i.e. alge- 
braic, correlations. This includes, as one particular example, 
the random field disorder, which will be discussed separately. 

We now give the results for these fixed points, first for 
short-range elasticity, then for LR elasticity, and compare with 
available numerical and exact results. The most important 
quantity to compute is the roughness exponent Since we 
have shown that X in j3.43> is universal to dominant order 
this proves universality of ( to the order in e studied here 
(i.e. 0(e^)). For LR disorder and for periodic fixed points 
we can also compute the universal amplitudes for the correla- 
tion function of displacements, and discuss their dependence 
on large scale boundary conditions. Anticipating a bit, let us 
summarize the general result that we use in that case, which is 
derived in Section IVll The T = disorder-averaged 2-point 
function for q ^ 0, q/m fixed, reads for any dimension d, in 
Fourier 

= {2T:f 6''{q + q')C{q) (4.1) 
C{q) = C{q - 0)Frf(q/m) (4.2) 
C{q = 0) = l{d)m-'^~'^'^ (4.3) 

The amplitude c{d) is given by the relation (exact to all orders 



in the present scheme): 

~c{d) 



(4.4) 



It is found to be universal only for long range and periodic 
disorder The scaling function, computed in Section fvTl for 
SR and LR elasticity, is always universal (independent of short 
scale details) and satisfies Fd{0) ~ 1 and 



Fd{z) ^ Bz-'^'^+^'^'l forz 
B = l + be + 0{e^) 



(4.5) 
(4.6) 



where b is computed in Section lVTI This gives us all we need 
for a calculation to O(e^) of the universal amplitude, e.g for 
the propagator in the massless limit m ^ q: 



C'iq) 
c{d) 



c{d){l + be + 0{^)) 



(4.7) 
(4.8) 



The result for C{q ~ 0) in presence of a mass is also in- 
teresting since it gives the fluctuations of the center of mass 
coordinate for an interface physically confined in a quadratic 
well. Although that situation would be interesting to study nu- 
merically, most numerical results are for finite size systems of 
volume L'^ (and m — + 0). We thus also define in that case: 



CL{q) = c'id)q-^''+'i^gdiqL) 



(4.9) 



with lim^^oo .^dl^) = 1. For peno(i/c boundary conditions 
q ~ 2TTn/L, n E Z"^ and n ^ 0. The prime indicates that the 
value of this amplitude depends on the large scale boundary 
conditions, i.e. it depends on whether e.g. a mass is used or 
periodic boundary conditions as an infrared cutoff. The ratio, 
computed in Section lVTl for short range elasticity. 



c{d) 



= 1 - 1.46935C+0(e2) 



(4.10) 



is unity only for periodic disorder, in which case the amplitude 
is independent of both large and small scale details. 

Before studying the different fixed points, let us mention 
an important property, valid under all conditions: If R{u) is 
solution of (I3.43> . then 



R{u) := k'^R{u/k) 



(4.11) 



is also a solution (for k a constant independent of m). We 
can use this property to fix R{0) or i?"(0) in the case of 
non-periodic disorder. (For periodic disorder the solution is 
unique, since the period is fixed.) 

A. Non-periodic systems: Random bond disorder 

Let us now look for a solution of our 2-loop FRG equation 
which decays exponentially fast at infinity as expected for SR 
random-bond disorder To this aim, we have to solve order by 
order in e the fixed-point equation ( I3.43> numerically. Making 
the ansatz 



R{u) = eri(ii) +6^2(1*) + 
C = eCi + e'C2 + . . . , 



(4.12) 
(4.13) 
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Taylor-expansion up to order 25 below: 



the partial differential equation to be solved at leading order 
is 



= (1 - 4Ci)ri(«) + Ciur[{u) + -r'Huf - r;'(w)r'/(0) 



1 = ri(0) 



(4.14) 



where we have used our freedom to normalize i?(0) :— e. 
i4.14\ has a solution for any , but only for one specific value 
of (i does this solution decay exponentially fast to 0, without 
crossing the axis, see figure^] The strategy is thus the fol- 
lowing: One guesses (i, and then integrates from to 
infinity. In practice, however, there are numerical problems 
for small u. One strategy, which we have adopted here, and 
which works very well, is to use the value of (i, to gener- 
ate a Taylor-expansion about 0. This Taylor-expansion is then 
evaluated at 0.5, where the numerical integration of (I4.14> is 
started, both forwards to infinity (which in practice is cho- 
sen to be 25) and backwards to 0. This enables to control 
the accuracy of both the Taylor-expansion and the numerical 
integration. The result for the best value 



ri{u) 



1-0.288797 7/2+0.0967487 ^^-0.0109959 
+0.000197282 7/^+0.0000162077 
+1.37054 10"^ u^+1.06127 10"^ 
+5.84538 10"^ u^-1. 50021 10" 
1.19821 10"" 
3.9358410- 
4.4915410" 
+6.77579 10" 
+4.1934810" 
+7.7804410" 
-4.3755710 
-6.74331 10 



-10 


u"-2.52931 10" 


11 


7/12 


-12 


1*13-4.9071710" 


13 


7/" 


-14 


1*1^-1.2175810- 


15 


7/16 


-16 


m"+2. 11465 10" 


16 


7/l« 


-17 


1*1^+6.49482 10" 


18 


7.2" 


-19 


1*21+5.52691 10" 


20 


7.22 


-21 


1*23-2.72231 10" 


21 


7/24 


-22 


7/25 + 0(^26) 







At second order in e, we have to solve 



(4.16) 



r2{u) - 4C2ri(7*) - 4Cir2(7i) + uC2r[{u) + 7tCir2(7*) 
+ 1 {r'iiu) - r'lm r'I'iuf - iri'(7/)rl"(0+)2 (4.17) 







2 

^2(0) 



(4.18) 



Ci = 0.20829806(3) (4.15) 



is given on figure [lO] (Note that in i46ll only the first four 
digits were given.) On this scale, Taylor-expansion and nu- 
merical integration are indistinguishable. The error-estimate 
on the last digit comes from moving the starting-point of the 
numerical integration (which was 0.5 above) up to 1, which 
allows for a crude estimate of the error. We also reproduce the 



where the last equation reflects our choice of R{0) = e. Note 
that to solve the 2-loop order equation, one has to feed in the 
solution at 1-loop order, both the Taylor-expansion about 
and the numerically obtained solution for larger u. Again (2 
is determined from the condition that the solution decays at 
infinity. Following the same procedure as at 1-loop order, we 
find 

C2 = 0.006858(1) . (4.19) 
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The function r2 is plotted on figurefH] The Taylor-expansion 
up to order 25 about reads 

r2{u) ^ -0.0604942 1*2+0.0345276 14^-0.00628098 
+0.000239628 w'V0.000019823 
+1.42202 10^^ u^+5. 17941 10~* 
-8.64456 10"^ u^-2. 72755 10~^ u^" 
^"m"-6. 23531 10 



-4.7860710 

-5.49541 10"'^ u"-8. 78473 10 
+1.30232 10~^^ u^^+3.60568 10 
+6.7239 10"^^ -u^+g. 51299 10 
+9.06111 10"" 1*1^-9.06201 10 



-1.5392210" 
-4.42973 10" 



^19 ^23. 



-21 ^25 



2.3656910 



-11 








-15 








-14 


u'^ 






-16 _ 


18 


U 


-20 








-19 






-20 







(4.20) 



One observes that Csr is necessarily bounded from above by 
e/4 as no SR solution can cross this value (to any order) with- 
out exploding. This reflect the exact bound for SR disorder 
6 < d/2, which simply means that optimization of energy 
must lower energy fluctuations compared to a simple sum of 
random numbers. Equality is obtained for the trivial constant 
eigenmode R{u) = i?.(0) corresponding to C = e/4, associ- 
ated to the fluctuation of the zero mode of the random poten- 
tial. 

We can now discuss our results for the roughness exponent. 
These are summarized in Table and compared to numer- 
ical simulations in d = 3,2 and the exact result for the di- 
rected polymer in = 1. A first observation is that the cor- 
rections compared to the 1-loop result have the correct sign 
and, further, that they improve the precision of the 1 -loop re- 
sult. Given the difficulties associated with this theory, this 
is a significant achievement. Second, the error bars given in 
Table [ni are estimated as half the 2-loop contribution, which 
should not be taken too literally, as it is difficult to obtain a 
good precision from only two terms of the series and no cur- 
rently available information about the large order behavior of 
this novel e-expansion. Third, one may try to improve the pre- 
cision using the exact result ( = 2/3 in = 1. Estimating the 
third order correction in the three possible Fade's in order to 
match — 2/3 for e = 3, we obtain consistently the values 
quoted in the fourth column of Table 1 121 We hope that these 
predictions can be tested in higher precision numerics soon. 

B. Non-periodic systems: random field disorder 

Let us first recall that at the level of the bare model the static 
random field disorder correlator obeys ~ — o'Im] at large 
|u| |46, 70], where a — (e/i)cr is proportional to the ampli- 
tude of the random field. 

If one studies the large u behavior in the ERG equation 
( I3.43> one clearly sees that the non-linear terms do not con- 
tribute, thus one has: 



Thus for a RE fixed point to exist, the O(e^) correction to ( 
has to vanish. 



Crf = e/3 . 



(4.22) 



This will presumably hold to all orders. Indeed it is clear that 
if there is a similar /9-function to any order, since each R car- 
ries at least two derivatives and at least one must be evalu- 
ated at u 7^ 0, the sum of all non-linear terms to a given 
finite order decreases at least as R" (u) ~ 1 /u. (This does 
not strictly excludes that summing up all orders may yield a 
slower decay, although it appears far fetched and does not oc- 
cur in the non-perturbative large-A^ limit.) The above value 
of C ensures that nfR{u) ^ — cr|u| in the effective action, i.e. 
non-renormalization of <j. 

Note that this argument based on long-range large u be- 
havior is a priori valid for any A. Since it is made on the R 
equation (no such argument can be made on the equation for 
A) it uses the property of potentiality. However, from ( I3.46t 
with C = e/3 one sees that A 7^ 1 is incompatible with the ex- 
istence of a fixed point, even of a fixed point with a supercusp. 
Thus, the only way to satisfy potentiality for the static random 
field problem seems to have a unrenormalized, = e/3 and 
A = 1 (the previous discussion of potentiality in Section III.D 
assumed short-range disorder). 

This must be contrasted with the theory of depinning, where 
we found that: 



dcp 



1'^ 



0.14331e) 



(4.23) 



following from Adop = — 1 in ( I3.47> . Since in that case the 
RG-flow is non-potential, it is clear that no similar argument 
as above exists to protect the value — e/i. (The force cor- 
relator is short range). The conjecture of I57i] thus appears 
rather unphysical in that respect. 

1. Fixed-point function 

We first study the fixed-point equation for 

K{u) = -R"{u) = ^y{u) (4.24) 
y(0) = 1 (4.25) 

and later use the rescaling freedom to tune the solution to the 
correct value of a at large scale u. 

The 2-loop ERG equation (|^^} becomes (A ~ 1): 



1 



2\ll 



(4.26) 



\{y'\y-i))" -\y'{o^fy" 



One can then integrate once with respect to u\ 

= uy-y'{y- 1) 



(4.27) 



\{y'\y~i))'~\y'{o^?y' 



-mdri 



(e-3C)a 



(4.21) 



There is no integration constant here because the second line 
precisely vanishes at it = 0+ (absence of supercusp). 
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one loop 


two loop 


estimate 


improved estimate 


simulation and exact 


d = 3 


0.208 


0.215 


0.215 ± 0.004 


0.214 


0.22 ±0.01 r82J 


d = 2 


0.417 


0.444 


0.444 ± 0.015 


0.438 


0.41 ±0.01 ^821 


d= 1 


0.625 


0.687 


0.687 ±0.03 


2/3 


2/3 



FIG. 12: First column: Exponents obtained by setting e = 4 — d in the 1-loop result. Second column: Exponents obtained by setting e = 4 — d 
in the 2-loop result. Third column: errors bars are estimated as half the 2-loop contribution. Fourth column: Improved estimates using the 
exact result ^cq = 2/3 in d = 1 (see text). 



The 1-loop solution involves the first line only. Dividing by 
y and integrating over u yields: 



yi - 1 - In yi 



(4.28) 



i.e. an implicit equation for y, which defines y — yi(u). It 
satisfies 

yi(0) = l, y[{0+) = -l 

y'l{o+) = I , vTio^) = -I ■ (4.29) 

We can put the 2-loop solution under a similar form. Making 
the ansatz 



— ^y~l~lny--F{y), (4.30) 



one obtains 

F{y{u))^l r^{y'\y-l)-y)' . (4.31) 

^ JO y 

At this order, one can replace y by j/i, i.e. use uy = y'{y — 1) 
to eliminate y' . This gives, changing variables from u to y: 



2 7i y dy V y-l 



The last term in the brackets is easily integrated. For the re- 
maining terms, we integrate by part and use ( 14.301 1 to replace 
by ?/ — 1 — In y: 

y-l-lny y(y-\-\ny) 1 _ 
F{y)^ Ay + = — T^lny. 

(4.33) 

This yields the final result 

F{y) = 2y - 1 ± ^ - \ Iny ± Li2(l - y) (4.34) 
1 - y 2 



(4.35) 



We find: 



F{V) -§(?/- 1)' - 3^(2/ - 1)^^ + 0((y - \f) (4.36) 

has a quadratic behavior around y = 1, similar to the 1-loop 
result, and corrects the value of the cusp. 



2. Universal amplitude 

Since we know the exact fixed point function up to a scale 
factor, we can now fix the scale by fitting the exact large |u| 
behavior to R(u) ^ —<^\ij\ where cr is the amplitude of the 
random field. The general fixed point solution reads: 



where ^ can be related to a as: 

/>oo 

o — I Au A(u) 

^0 

We need 



71 



eh- 



Auy(u) = / dyu(y) 
/o Jo 

= 71 + £72 

/ dy V2(y-i-lny) 
F{y) 



(4.37) 



(4.38) 



= 0.775304245188 

12 = - I dy 



JO 3^2(7 
-0.13945524. 



(4.32) One can now express 



3(7 



l-lny) 



2/3 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



and thus compute, using ( 14.41 1 the universal amplitude ( 14.31 1 
associated to the mode q — Q'm presence of a confining mass: 



_L 

%d) = (t3 ' (71 ± £72)^^(6/1)^3 



(4.43) 



(0 ' (71 + £72)" 



_(47r)d/2 



where one has restored the factors e/i absorbed in A and a. 
Expanding all factors in a series of e one finds: 

c{d) = ei/3(3_52459 o.725079e + ©(e^))^^/^ , (4.44) 



The lowest order was obtained in Ref. 17011 and we have ob- 
tained here the next order corrections. It is interesting to com- 
pare our result with the exact result in d 0, which is I69H : 



c[d = 0) = 1.05423856519 . . .cr^/^ 



(4.45) 
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While the simple extrapolation setting e = 4 of ( I4.44> to one 
loop c{d = 0) = 5.59(7^/'^ is very far off, to two loop it gives 
c{d = 0) = 0.99cr^/^, surprisingly close to the exact result. It 
was noted in Ref. |70] that extrapolation of the 1-loop result 
could be considerably improved by not expanding ( I4.43t in 
e but instead directly setting e = 4 (with 72 = 0) in ( I4.43t . 
That gives ci{d = 0) = 0.821(t^/'', an underestimate already 
reasonably close from the exact result. We extend this pro- 

— 2/3 

cedure to two loop by truncating the e expansion of ly 
to second order in (I4.43> . and then set e = 4. This yields 
C2{d = 0) = 1.22(7^/^, and the exact result is then halfway 
between ci{d — 0) and C2{d = 0). To summarize, our 2-loop 
corrections J4.44> have the correct sign and order of magni- 
tude to improve the agreement with the exact result in d = 0. 

The universal amplitude for the massless case ( 14. 7> (or q ^ 
to) is obtained from ( 14. 8> with b = —1/3 from Section lvTl as: 



1 



c{d) = 5(d) 

= e^/^ [3.52459 - 1.89994e + 0{e^)\ a^'^ , (4.46) 



and writing c{d) — c{d)/{l + ^e) should provide a reason- 
able extrapolation to low dimensions. Finally, we recall that 
for random field disorder, this coefficient is different for dif- 
ferent large scale boundary conditions. The result for periodic 
boundary conditions can be obtained from formula ( I4.10t . 

In Ref. fTd\, the 1-loop result was compared to the result of 
the Gaussian Variational Method (GVM). It is instructive to 
pursue this comparison to two loops. We get from ItoII : 



CGVM(d) 

cgvm(c^) 

CGVM {d) 



1 




y 


_(4^)d/2_ 



1 



-(ja 



(4.47) 



= .1/3 



ei/-'(3.69054 - 0.894223e + 0{e^))a 



2/3 



2C 



1 



2 



2C\ ^(e-2C)/2 



sin(7r(e-2C)/2) 
(4.48) 



where in the last line we have inserted C = e/3 and performed 
the e expansion. Thus one finds, quite generally that b^ar = 
36/4. As noted in |701 to one loop the FRG and the GVM 
give rather close amplitudes (differing by about 5 per cent). 
We see here that to two loop, i.e. next order in e, the difference 
increases. Finally, 

CGVM(d) = e^/^(3. 69054 - 1.81686e + 0(e^))cr^/^ 
and the coefficient remains rather close to the one in ( I4.46I I. 

C. Generic long range fixed points 

There is a family of fixed points such that 

R{u) - (4.49) 

associated with 



c 



2(1 + 7) 



(4.50) 




FIG. 13: The fixed-point function yi(u) at 1-loop order for non- 
periodic disorder. 



These fixed points where found for infinite N in any d in 
Ref. i36[ It^ (we use the same notations). They were stud- 
ied to first order in e for any N in f47], and argued to be stable 
only for 7 < 7*((i) the value of the crossover to short range 
identified in |47] as CsR = C,i.k{i* {d)) . 

Here, we have not studied these fixed points in detail but we 
note that the 2-loop corrections do not change C,, by the same 
discussion as for the random field case 7 = 1/2. They will 
however affect the amplitudes. 

D. Periodic systems 

1 . Fixed point function 

For periodic R{u) as e.g. CDW there is another fixed point 
of ( I3.43t . It is sufficient to study the case where the period 
is set to unity, all other cases are easily obtained using the 
reparametrization invariance of equation ( 14.1 1> . No rescaling 
is possible in that direction, and thus the roughness exponent 
is 



C = o 



(4.51) 



The fixed-point function is then periodic, and can in the in- 
terval [0, 1] be expanded in a Taylor-series in u(l — u). Even 
more, the ansatz 

k{u) = [aie + 026^ + ...) + (6ie + 626^ + • ■ • ) u^(l - uf 

(4.52) 

allows to satisfy the fixed-point equation ( I3.43t to order and 
will presumably work to all orders. For a more general case 
of this see Ref. |68]. 

To gain insight into the more general case, let us write the 
fixed point for ( I3.43> with arbitrary A: 



R*{u) = — ^ + (3 - 2A)-^ 
^ ' 2592 ^ ' 7776 



+ (A- 1) uil 

^ ^^432 ^ 



u) 



W2 108 / ^ ' 



(4.53) 
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One can see on this solution that A = 1 is the only value 
which avoids the appearance at two loops of the supercusp, 
i.e. a cusp in the potential correlator R{u) rather than in the 
force correlator A(ii). 

The same discussion can be made on the the flow equation 
of A(u) by taking two derivatives of ( I3.43> . One finds that 
there is a priori an unstable direction corresponding to a uni- 
form shift in A(it) A(ii) + est. While this is natural in e.g. 
depinning, it is here forbidden by the potential nature of the 
problem which requires 



duA{u) = 



(4.54) 



since in a potential environment, the integral of the force over 
one period must vanish. This is indeed satisfied for the fixed 
point for A{u) 

A*{u) = -R*"{u) 

(4.55) 



only if A = 1: 



1 



„ d.A*(.).-(A-l) 



(4.56) 



The values for depinning are obtained by setting A = —1: in 
that case, the problem becomes non-potential at large scales. 

2. Universal amplitude 

This fixed point impUes for the amplitude of the zero mode in 
presence of an harmonic well, defined in ( 14. 3> . using ( 14. 4> : 

= 2.19325e - 0.680427e^ + 0{e^) (4.57) 

In the other limit m <^ q one obtains the amplitude (using 
b = -1 from SectioniyH: 



c(d) 



er(|) 1^36 ^ 108 
2.19325e - 2.87367e2 + 0{e^) 



(4.58) 
(4.59) 



Note that we prove in Section ( IVU that this amplitude is inde- 
pendent of large scale boundary conditions, and is thus iden- 
tical for e.g. periodic boundary conditions and in presence of 
a mass. As can be seen from ( I4.10t this is a consequence of C, 
being zero. 

This can be compared to the GVM method j2^|3- 



CGVM(d) = (4 - d)2'''\^-^r ^-j (4.60) 
= 2e- 2.9538e2 + 0(e3) (4.61) 



with coefficients surprisingly close to the e-expansion. 

It is interesting to compare predictions in d = 3. We recall 
that we are studying a problem where the period is unity, the 
general case being obtained by a trivial rescaling in u. Since 
( I4.59> has a poor behavior (and so does ( 14. 6H which resums 
into ( I4.60> '). it is better to use instead ( I4.58> . It was indeed 
noted in |26^,'23l that the improved 1-loop prediction ci{d ~ 
3) obtained by setting e = 1 and ignoring the e^/108 term 
in ( I4.58> yields a value rather close to the prediction of the 
GVM: 

ci(d = 3) = 27r/9 = 0.6981 (4.62) 
CGVM(d = 3) = 1/2. (4.63) 

Including the 2-loop e^/108 term now gives C2{d = 3) = 
0.4654 and C2{d = 3) = 0.5235 for the two Fades respec- 
tively. This type of extrapolation makes the GVM and FRG 
predictions get closer when including the 2-loop corrections. 
On the other hand, comparison of ( I4.59> and ( 14. 6U suggests 
thatc(o?) > CGVM(rf)- 

This is in reasonable agreement with the numerical results 
of Middleton et al. |84|. They obtained good evidence for 
the existence of the Bragg glass (i.e. its stability with respect 
to topological defects predicted in |26, 2^). They measure 
directly the correlation ( 14. 7> and obtain strong evidence for 
the behavior ( I4.58t (as well as the correct correction to scaling 
behavior) with 



2c(d = 3) wl.04, 



(4.64) 



(their amplitude A is twice our c{d)) which lies in between 
the GVM and the 1-loop FRG. (More precisely two different 
discretizations gave 2c{d = 3) = 1.01 ± 0.04 and 2c{d = 
3) = 1.08 ±0.05). 

Another interesting observable is the slow growth of dis- 
placements characteristic of the Bragg glass: 



(4.65) 



at large x. Performing the momentum integral from i4.7i . one 
obtains: 



A, 



{iTT)d/^r{d/2) 

4sin(7re/2) 



c{d) 



7re(l-|) [36 108 
If one expands each factor in e it yields: 



A„ 



18 ^ 108 



For comparison, the GVM gives 



(4.66) 



(4.67) 



(4.68) 



Here extrapolation directly setting e = 1 in (I4.67> looks pos- 
sible, and yields A3 — 0.0556 to one loop increasing to 
A3 = 0.0648 to two loop. On the other hand, setting e — 1 
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in (I4.66> yields instead ^3 = 0.0707 to one loop decreas- 
ing to = 0.047 at two loops. The GVM gives the result 
^3,GVM = 0.0507. 

Another interesting observable is: 



w"^ = Bd In L 



w 



u~ 



1 

1? 



(4.69) 
(4.70) 



where L is the linear system size. In Ref. 18411 it was assumed 
that Bd = Ad/2 thus in d = 3, B3 = c(3)/(27r2) yielding 
a value of c(d) consistent with the direct measurement of this 
quantity 1 104|. This was also done in |85] where it was de- 
duced from a measurement of Bd that 0.98 < 2c{d = 3) < 
1.11 11 05.1 . Although this is a reasonable approximation, it 
is not exact. Indeed the quantity Bd, contrary to c(d), de- 
pends on the (large scale) boundary conditions. It is of course 
universal, since it does not depend on small scale details. Its 
value can be computed e.g. for periodic boundary conditions 
and pinned zero mode, and depends on the whole finite size 
scaling function ( 14. 9> computed in Section lvTI 



-(d+2C 



9d{qL) 



(4.71) 



As shown recently, nP' fluctuates from sample to sample and 
the full distribution P{w'^) averaged over disorder realizations 
was computed for the depinning problem LS4,,,S2ll- 

E. Long range elasticity 

Let us now consider the case of long range elasticity. There 
are physical systems where the elastic energy does not scale 
with the square of the wave-vector q as -Eciastic ~ but as 
inelastic ~ l^l". In this situation, the upper critical dimension 
is dr = 2a and we define: 



■.= 2a-d. 



(4.72) 



The most interesting case, a priori relevant to model a contact 
line is a = 1, thus dc = 2. For calculational convenience, we 
choose the elastic energy to be 

E^ustic - {q^ + m^)^ . (4.73) 

This changes the free correlation to: 

Gab{q) = Sab- 



' {q^ + m?) 2 
The energy exponent in that case is: 

e» = a - d + 2C . 



(4.74) 



(4.75) 



The changes are very similar to the case of Ref. \61V so we 
summarize them here only briefly. The /3-function is still 
given by but with the integrals replaced by: 



^1 ~ 



1 



(g2 



-,2\a 



-er(e/2) 
= m " \' / I e 



1 



(4.76) 



(4.77) 



and thus the /3-function is given by (I3.43> with: 

2 6(2/<°)-(/j"))2) 



(e/j"))2 
1 + it 

~" 

0(6). 



(4.78) 



(See appendix F of Ref. |y67il ). And of course the relation 
( !3.42> between R and R is identical except that e/i must be re- 
placed by el[°'\ Since X^"^ is finite, the /3-function is finite; 
this is of course necessary for the theory to be renormalizable. 
For the cases of interest a — 1 and a = 2, we find 



Xi^) = 1 
X^^^ = 41n2 . 



(4.79) 
(4.80) 



The exponent (as a function of e) and the fixed point function 
is thus changed only at two loops. 

Let us now give the results in the cases of interest: 

1. Random bond disorder 



The solution of ( I3.43> with X 
second order in e as: 



X^"' can be written, to 



R{u) = eri{u) + e^X^°'^r2iu) + . . . (4.81) 
C = eCi+e2x(")C2 + ... , (4.82) 

since the equation ( I4.17t for r2 (u) is linear Thus one has for 
any a: 

C = 0.20829806(3)e + 0.006858(l)X(")e2_^O(e3) . (4.83) 
For the case of most interest a — 1, X^^^ = 4 In 2 one finds: 



C 0.20829806(3)e+ 0.0190114(3)e2 
e = 2-d 



(4.84) 



and = 2C. 

It would thus be interesting to perform numerical simu- 
lations in d = 1 for the directed polymer with LR elastic- 
ity. This would be another non trivial test of the 2-loop cor- 
rections. The 1-loop prediction is = 0.208, significantly 
smaller than the roughness for SR elasticity C — 2/3. The 
naive 2-loop result is (setting e = 1), C w 0.227 ± 0.01. Er- 
ror bars are estimated by half the difference between the 1- 
loop and 2-loop results. Note that the bound 9 < d/2 implies 
C < 1/4 in d = 1, already rather close to the 2-loop result. 

2. Random field disorder 

The exponent is still 



c 



(4.85) 



and was indeed measured in experiments on an equilibrium 
contact line |30]. It would be of interest to measure the uni- 
versal distributions there, such as the one defined in L86..87.1 . 
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The fixed point function is given by ( I4.3Q> and J4.34t upon 
replacing F{y) X^°'^F{y). The amplitude of the zero 
mode in a well c{d) is now given by: 



c{d) = a 



- ^2/3 



e\i/3 / , , \-2/3 



(J(a))^l/3 

(4.86) 
(4.87) 



and the amplitude of the massless propagator 

c{d) =S(d)(l + 6„e) . 
where ha is given in ( I6.14t setting = 1/3. 

3. Periodic disorder 

The fixed point becomes: 

A*(u) = -W"{u) 

^^ + |^x(")-(^4x(")).(l-.).(4.88) 

For the periodic case, the universal amplitude reads: 



and 

c{d) =5(d)(l + 6„e) . 
Setting Ci = in ( I6.14> yields 

4 



-dd) 



(47r)'i/2r(d/2) 
Using e — 2oL — d, this gives 

1 4X(") +3(7 + 7/'(q;)) +66„ 
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which in the case of a = 1 takes the simple form 
18 108 



(4.89) 



(4.90) 



(4.91) 



(4.92) 



(4.93) 



V. LIFTING AMBIGUITIES IN 
NON-ANALYTIC THEORY 

A. Summary of possible methods 

As we have seen above ambiguities arise in computing the 
effective action at the level of 2-loop diagrams if one uses a 
non-analytic action. One can see that these arise even at the 1- 
loop level for correlations (see below Section lVH . To resolve 
this issue, our strategy has been to use physics as a guide and 
require the theory to be renormalizable, potential and without 
supercusp. This pointed to a specific assignment of values to 
the "anomalous" graphs. The physical properties of the en- 
suing theory, studied in the previous section, were found to 
be quite reasonable. Of course, one would like to have a bet- 
ter, more detailed justification of the used "prescription". Al- 
though we do not know at present of a derivation of this theory 



from first principles, we have developed a set of observations 
and a number of rather natural and compelling "rules" which 
all lead to the same theory. We describe below our success- 
ful efforts in that direction as well as some unsuccessful ones, 
which illustrate the difficulty of the problem. 

A number of approaches can be explored to lift the ambi- 
guities in the non-analytic theory. We here give a list; some of 
the methods will be detailed in the forthcoming sections. 

1) Non-zero temperature: At T > previous Wilson 1- 
loop FRG analysis fss'.'s?,"?^,'!?] found that the effective ac- 
tion remains analytic in a boundary layer u ^ T. However, 
since the rescaled temperature ( I2.24> flows to zero as T ~ 
as m ^ (temperature being formally irrelevant) all (even) 
derivatives of R{u) higher than second grow unboundedly as 
TO ^ 0, for instance i?""(0) - R*"'{{)+flf (in terms of the 
zero temperature fixed point function). On a qualitative level 
one can thus see how finite T diagrams such as E in Fig.|3l 
yielding 



TR""{0)R"{u) -> R*"'{0+)^R"{u) 



(5.1) 



can build up "anomalous" terms in the /3-function, hence con- 
firming what is found here 1 70]. However, correctly and quan- 
titatively accounting for higher loops is a non-trivial problem 
as stronger blow-up in 1 /T'' seem to arise. In fact each new 
loop brings two derivatives and a propagator, hence an addi- 
tional factor 1/T. Despite some recent progress, a quantitative 
finite-temperature approach which would reproduce and jus- 
tify the present e expansion has proved difficult|74, 75]. Not 
only for technical reasons, as methods using exact RG where 
found to be appropriate, but also for physical reasons, as an 
extension to non-zero T must also handle low-lying thermal 
excitations in the system (e.g. droplets). A theory from first 
principles at T > is thus presently not available and will 
not be further addressed here. All other methods use a non- 
analytic action. 

2) Exact RG: Exact RG methods directly at T = have 
been studied to one loop ]70, 89] and two loopsf7l','9^. Al- 
though it does yield interesting insights into the way to handle 
ambiguities (see below), and confirm the present results, it 
suffers from basically the same problems as described here. 

3) Direct evaluation of non-analytic averages: In this ap- 
proach one attempts a direct evaluation of non-analytic av- 
erages (e.g. in fully saturated diagrams). For instance, ex- 
panding at each vertex the disorder R{u^ — u^) in powers of 
l^a ^ I using the proper non-analytic Taylor expansion: 



R"{u) = i?"(0) + i?"'(0+)|u| + R""{0+)u^ + 



(5.2) 



one can try to compute directly all averages in vertex functions 
and correlations. After performing a few Wick contractions 
one typically ends up with averages involving sign functions 
or delta functions. These can be computed in principle using 
the free Gaussian measure. For instance, using formulae such 
as: 



(sgn(u)sgn(i;)). 



(5.3) 
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Although promising at first sight, the results are disappoint- 
ing. Averages over the thermal measure involve many changes 
of signs which kill all interesting divergences indicating that 
some physics is missing. The method, briefly described in Ap- 
pendix|B]is thus not developed further. A dynamical version 
of this method which is similar in spirit i66l l67ll . did work 
for depinning, although there it simply identified with another 
method used below, the background field (which, for depin- 
ning is Uxt ^ vt + Uxt see below). 

4) Calculation of T{u) with excluded vertices and sym- 
metrization: A valid, general and useful observation (not lim- 
ited to this method) is that if one uses the excluded vertex 

(5.4) 

then all Wick contractions can be performed without ambigu- 
ities. The excluded vertex is as good as the non-excluded one 
since one can always add a constant — ni?(0) to the action of 
the model \2.2\ . Thus one can compute without any ambiguity 
the effective action T{u) for an "off-diagonal" field configu- 
ration 

M° such that ^ for ah a ^ 6 , (5.5) 

since then no vertex is ever evaluated at u = 0. The drawback 
is that one ends up with expressions containing terms such as 

^ R"{Ua--Ub)R"'{Ua-Ub)R"'iub-Uc) , (5.6) 

which superficially looks like a three replica term, but due to 
the exclusions, may in fact contain a 2-replica part which can 
in principle be recovered from the above by adding appro- 
priate diagonal terms, using that p-replica parts are properly 
defined as free replica sums e.g. from a cumulant expansion. 
The 2-replica part of i5.6\ thus naively is 

^^i?"(o)i?"'(o)i?"'K-«e) 

ac 

+ ^i?"(Ua - Ub)R'"iUa - M6)i?"'(0) (5.7) 

ab 

and one is again faced with the problem of assigning a value 
to i?"'(0). The calculation with excluded vertices thus yields 
a sum of p-replica terms with p > 2 and to project them onto 
the needed 2-replica part, one may need to continue these ex- 
pressions to coinciding arguments u° = i^ . 

The symmetrization method attempts to do that in the most 
"natural" way. Using the permutation symmetry over replicas 
and the hypothesis of no supercusp yields a rather systematic 
method of continuation. Surprisingly, it fails to yield a renor- 
malizable theory at two loops. We identified some difference 
with methods which do work, but the precise reason for the 
failure in terms of continuity properties remains unclear. It 
may thus be that there is a way to make this method work but 
we have not found it. Being interesting in spirit this method is 
reported in some details in AppendixIXI 



If one renounces to the projection onto 2-replica terms one 
can, in a certain sense, obtain renormalizability properties. 
This generates an infinite number of different replica sums and 
seems to be not promising, too. It is described in Appendix|F| 

We now come to methods which were found to work, and 
which will be described in detail in the next section. In all 
of them one performs the Wick contractions in some given 
order (the order hopefully does not matter) and uses at each 
stage some properties. The fact that one can order the Wick 
contractions stems from the identity, which we recall, for any 
set of mutually correlated Gaussian variables Ui. 

{u^W{u)) = ("'"j) (^«.^(^)) (5-8) 

under very little analyticity assumption for W{u\ which can 
even be a distribution. At each stage one can either use ex- 
cluded or non-excluded vertices as is found more convenient. 

5) Elimination of sloops: We found another method, 
which seems rather compelling, to determine the 2-replica 
part of terms such as ( I5.6l l. It starts, as the previous one, 
by computing (unambiguously) diagrams with the excluded 
vertices. Then instead of symmetrization, one uses identities 
derived from the fact that diagrams with free replica sums and 
which contain sloops cannot appear in a T = theory and 
can thus be set to zero. Further contracting such diagrams 
generates a set of identities which, remarkably, is sufficient 
to obtain unambiguously the 2-replica projection without any 
further assumption. It works very nicely and produces a renor- 
malizable theory, as we have checked up to three loops. In 
some sense, it uses in a non-trivial way the constraint that we 
are working with a true T = theory. This method is detailed 
below. 

6) Background field method: This method is similar to 
method number 3 except that the vertex R(u) at point x is 
evaluated at the field = + then expanded in 
which then are contracted in some order. This amounts to 
compute the effective action in presence of a uniform back- 
ground field which satisfies j5.5> . Thanks to this uniform 
background and upon some rather weak assumptions, the am- 
biguities seem to disappear. The method is explained below. 

7) Recursive construction: An efficient method is to con- 
struct diagrams recursively. The idea is to identify in a first 
step parts of the diagram, which can be computed without am- 
biguity. This is in general the 1-loop chain-diagram ( 13. U . In 
a second step, one treats the already calculated sub-diagrams 
as effective vertices. In general, these vertices have the same 
analyticity properties, namely are derivable twice, and then 
haveacusp. (Compai-e with (i?"(u)-i?"(0))i?"'(u)2- 
i?"(it)i?"'(0+)'^). By construction, this method ensures renor- 
malizability, at least as long as there is only one possible path. 
However it is not more general than the demand of renormal- 
izability diagram by diagram, discussed below. 

8) Renormalizability diagram by diagram: In Section III 
we have used a global renormalizability requirement: The 1- 
loop repeated counter-term being non-ambiguous one could 
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fix all ambiguities of the divergent 2-loop corrections. How- 
ever, as will be discussed in |73], this global constraint ap- 
pears insufficient at three loops to fix all ambiguities. Fortu- 
nately, one notes that renormalizability even gives a stronger 
constraint, namely renormalizability diagram by diagram. 
The idea goes back to formal proofs of perturbative renormal- 
izability in field-theory, see e.g. (H ill HI H 11 IH IH 
^8]. These methods define a subtraction operator R. Graph- 
ically it can be constructed by drawing a box around each 
sub-divergence, which leads to a "forest" or "nest" of sub- 
diagrams (the counter-terms in the usual language), which 
have to be subtracted, rendering the diagram "finite". The 
advantage of this procedure is that it explicitly assigns all 
counter-terms to a given diagram, which finally yields a proof 
of perturbative renormalizability. If we demand that this proof 
goes through for the functional renormalization group, the 
counter-terms must necessarily have the same functional de- 
pendence on R{u) as the diagram itself. In general, the 
counter-terms are less ambiguous, and this procedure can thus 
be used to lift ambiguities in the calculation of the diagram 
itself. By construction this procedure is very similar to the 
recursive construction discussed under point 7. 

It has some limitations though. Indeed, if one applies this 
procedure to the 3-loop calculation, one finds that it renders 
unique all but one ambiguous diagram, namely 



(5.9) 



which has no subdivergence, thus there are no counter-terms, 
which could lift the ambiguities. Thus this diagram must be 
computed directly and we found that it can be obtained unam- 
biguously by the sloop elimination method|73|. 

9) Reparametrization invariance: From standard field 
theory, one knows that renormalization group functions are 
not unique, but depend on the renormalization scheme. Only 
critical exponents are unique. This is reflected in the freedom 
to reparametrize the coupling constant g according g — > g{g) 
where g{g) is a smooth function, which has to be invertible in 
the domain of validity of the RG /^-function. 

Here we have chosen a scheme, namely defining R{u) from 
the exact zero momentum effective action, using dimensional 
regularization, and a mass. One could explore the freedom 
in performing reparametrization. In the functional RG frame- 
work, reparametrizations are also functional, of the form 




R{u) 



R{u) = R[R]{u) 



(5.10) 



Of course the new function R{u) does not have the same 
meaning as R{u). Perturbatively this reads 

R{u) — > R{u) = R{u) + B{R,R){u) +0{R^) , (5.11) 

where B{R, R) is a functional of R. For consistency, one 
has to demand that B{R, R) has the same analyticity proper- 
ties as R, at least at the fixed point R = R*, i.e. B{R,R) 
should as R be twice differentiable and then have a cusp. 
A specifically useful candidate is the 1-loop counter-term 



can convince oneself, that by choos- 
ing the correct amplitude, one can eliminate all contributions 
of class A, in favor of contributions of class B. Details can be 
found in li73ll . 

Apart from methods 3 and 4 which did not work for reasons 
which remain to be better understood, methods 2,5,6,7,8,9 
were all found to give consistent result, making us confident 
that the resulting theory is sufficiently constrained by general 
arguments (such as renormalizability) to be uniquely identi- 
fied. Let us now turn to actual calculations using these meth- 
ods. 



B. Calculation using the sloop elimination method 

1. Unambiguous diagrammatics 

Let us redo the calculation of Section llll Bl using excluded ver- 
tices. From now on we use sometimes the short-hand nota- 
tions 



Rab = R{U^ ~ U^) 



R 



(p) 



(5.12) 



whenever confusion is not possible. 

The resulting diagrammatics looks very different from the 
usual unexcluded one. When making all four Wick contrac- 
tions of the 2-loop diagrams A, B and C in Fig. |6l between 
three unsplitted vertices one now excludes all diagrams with 
saturated vertices, but instead has to allow for more than 
two connected components and for sloops. The splitted ex- 
cluded diagrams corresponding to classes A, B and C are 
given in Fig. [H] There is an additional multiplicative coef- 
ficient l/(mi!TO2!'7i3!TO4!) in the combinatorics for each pair 
of unsplitted vertices (say ab and cd) linked by an internal 
line where mi propagators link ac, rn2 link ad, link he, 
1712 link hd. (This is equivalent to assigning a color to each 
propagator). 

Let us denote by 5T — —^^S'^^R the 2-loop contribution 
of all diagrams of class A to the effective action. One finds: 



R — r \^ D" f T}'" \2 , nil nlll nin 

R- [2^Rab[Rab) + 2^ Rab^abR-ac 



a^b 
1 



E 

a^b.a^c.b^c 



jD/f in/ff -nff/ 
2 / ^ ^ab^ac-^bc 



/// \ 2 

ac) 



a^b.a^c 



O- \^ P'" f?"' 

9 ^ab^ac^ad 



a^b .a^c.a^d 



Ia 



(5.13) 



coming respectively and in the same order from graphs a, j3, 
7, 5+rj (they are equal) and A in Fig. The only graph 
common to excluded and free-sum diagrammatics is a which 
is graph b of Fig.0 since all the other graphs in Fig. Q have 
saturated vertices. 
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FIG. 14: 2-loop diagrams corresponding to Larkin's hat (top) and 
banana (bottom). The graph a is a 1-replica term, /3, 7, 5 and rj are 
improper three rephca terms and A is an improper four replica term. 



Similarly, the graphs of class B give a total contribution: 



1 \ ^ p//// p" p" I ^ \ ^ p"" p" p" 

2 2^ ^ab "-ab^ab + ^ ^ab ^ac^a: 

a^b a^b,a^c^a^d 
1 jy//rr T>ii nil , nim p" p" 



a^b^a^c^b^d 



a^b^a^c 



(5.14) 



coming respectively and in the same order from graphs a', (3' , 
7', (5' in Fig. [21 Again, the only graph common to excluded 
and free-sum diagrammatics is a' which is graph g of Fig.|8] 
since all the other graphs in Fig.|8]have saturated vertices. 

The contribution 5^ i? of the diagrams of class C is given 
in Appendix Icl Note that adding a tadpole does not alter the 
structure of the summations in the excluded-replica formal- 
ism, since a tadpole can never identify indices on different ver- 
tices. This indicates that class C does not contain a 2-replica 
contribution, but starts with a 3-replica contribution (times T). 
This is explained in more details in AppendixIDl 

One can first check that when R{u) is analytic one recovers 
correctly the same result as (I3.27> setting the last (anomalous) 
term to zero. Adding and subtracting the excluded terms in 
(I5.13> to build free replica sums (using R"'{0) = in that 



case), or equivalently lifting all exclusions but replacing ev- 
erywhere 



R 



(p) 

ab 



(5.15) 



and then expanding and selecting the 2-replica part, one finds 
the contributions 



a R"{u)R"'{uf 
\R"{Q)R"'{uf 



7 



S + 1] 



-R"{0)R"'{uf 



(5.16) 



Similarly in ( 15.141 1 one obtains 



1 



^ _ tR""(u)R"{uf 

(3' ^-f' ^ ^R""{0)R"iO)R"{u) + ^R"{OfR""{u) 

5' ^ -R""{Q)R"{Q)R!'{u) - R"{Q)R""{u)R"{u) . 

(5.17) 

We now want to perform the same projection for a non- 
analytic R{u). 

2. The sloop elimination method 

The idea of the method is very simple. Let us consider the 1- 
loop functional diagram (a) in Fig.|2]which contains a sloop. 
It is a three replica term proportional to the temperature. In a 
T = theory such a diagram should not appear, so it can be 
identically set to zero: 



W 



— Y^R"{uf)R"{ul'^)^Q 

abc 



(5.18) 



It is multiplied by G{x — y)^, which we have not written. We 
will also omit global multiplicative numerical factors. Pro- 
jecting such terms to zero at any stage of further contractions 
is very natural in our present calculation (and also e.g. in the 
exact RG approach, where terms are constructed recursively 
and such forbidden terms must be projected out). It is valid 
only when (i) the summations over replicas are free (ii) the 
term inside the sum is non-ambiguous. These conditions are 
met for any diagram with sloops, provided the vertices have 
at most two derivatives. (One can in fact start from vertices 
which either have no derivative or exactly two.) 

Let us illustrate the procedure on an example. We want to 
contract W with a third vertex R at point z, i.e. we first write 
the product: 



de 



1 

rp4 



E 

a^b.a^CA 



R'libR'acRde = , (5.19) 



where implicitly here and in the following the vertices are at 
points X, y, z in that order. We will contract the third vertex 
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twice, once with the first and once with the second , i.e. look 
at the term proportional to G(a; — y)^G{x — z)G{y — z). Note 
that since we will contract each vertex, we are always allowed 
to introduce excluded sums (clearly the diagonal terms a ^ b, 
a ^ coi d = e give zero, since Rat and its two lowest deriva- 
tives at a = 6 are field independent constants). Performing the 
first contraction (i.e. inserting Sad — Sae — 5bd + 5be multiplied 
by the exclusion factors (1 — 5ab){^ — <5ac)(l — 5de) yields (up 
to a global factor of 2): 



1 

T3 



= 0. 



Ed'" nil pi _ nlll n" n' 

^ab^ac^ae / , ^ab^ac^be 

a^b.a^c.a^e a^b^a^c.b^e 

(5.20) 

Similarly, the second contraction then yields (up to a global 
factor of 4): 

— ri fl'" R" -I- \^ f?'" R'" R" 

\_2 / ^ -i^ab^ac^ae ' / , ^ab^ac^ac 



a^b,a^c^a^e 



a^b,a^c 



1 Tjlll Tjlll TD'I ^ TD'II T3'll p" 

^ab^ab^ae ^ ^ ^ab"-ac^bi 



a^b^b^e 



a^b^a^cM^c 



= . 
(5.21) 



This non-trivial identity tells us that the sum of aU terms (or di- 
agrams) generated upon contractions of diagram (a) of Fig.|2l 
(i.e. the 1-loop sloop-diagram equivalent to term W in ( I5.18» 
with other vertices, must vanish. Stated differently: A sloop, 
as well as the sum of all its descendents vanishes. Note that 
this is not true for each single term, but only for the sum. 

A property that we request from a proper p-replica term is 
that upon one self contraction it gives a (p — l)-replica term. 
It may also give T times a p-replica term (a sloop) but this 
is zero at T = 0, so we can continue to contract. Thus we 
have generated several non-trivial projection identities. The 
starting one is that the 2-replica part of ( I5.18> is zero, since 
(I5.18> is a proper 3-replica term. Thus, ( I5.19> prior to the 
exclusions, is a legitimate 5-replica term, and its 4-replica part 
is zero. Upon contracting once we obtain that the 3-replica 
part of ( I5.20> is zero. The final contraction tells us that the 
2-replica part of (I5.21> is zero. This is what is meant by the 
symbol "=" above and the last identity is the one we now use. 

Indeed compare ( 15. 2H with ( I5.13> . One notices that all 
terms apart from the first in ( I5.13> appear in ( I5.21> . and with 
the same relative coefficients, apart from the thkd one of 
(I5.13> . Thus one can use (I5.21> to simplify (I5.13> : 



5^'^R 



\^R'ab{R'abf 



E KbiKc? 



I A ■ (5.22) 



The function R"'{uY, which appears in the last term, is con- 
tinuous at u = 0. It is thus obvious how to rewrite this expres- 
sion using free summations and extract the 2-replica part 

5fR{u) = \{R"{u) ~ R"{0))R"'{uf 

-R"'{0+fR"{u)\lA, (5.23) 

which coincides with the contribution of diagrams A in (|^^} 
with A = 1. 



We can write diagrammatically the subtraction that has 
been performed 



sTr = 





(5.24) 



where the loop with the dashed line represents the sub- 
diagram with the sloop, i.e. the term ( I5.18> (with in fact the 
same global coefficient). The idea is of course that subtract- 
ing sloops is allowed since they formally vanish. 

There are other possible identities, which are descendants 
of other sloops. For instance a triangular sloop gives, by a 
similar calculation: 



a^b a^b.a^c 

. T}ii ( r>iii \2 . rtiii jyrri n 

+ Rbc(Rab) + 2^ Rac-n-bcR. 



a^b,b^c 



a^c^b^c^c^d 



(5.25) 



(2) 

This however does not prove useful to simplify S)^ R. 

Since the above method generates a large number of iden- 
tities, one can wonder whether they are all compatible. We 
have checked a large number of examples (see the 3-loop cal- 
culations in 1 73 1) and found no contradictions, although we 
have not attempted a general proof. 

The diagrams B and C are computed in Appendix IdI One 
finds by the same procedure 



SbR=\r""{u){R"{u) 
ScR = 0, 



R"iO)f 



(5.26) 
(5.27) 



confirming our earlier results in section llll B 2l and lIII B 31 

C. Background method 

In the background method, one computes r[u] to two loops 
for a uniform background u such that Uab 7^ for any a ^ b. 
We start from 



ipi 



(5.28) 



Taylor expand in v^, and contract all the v fields keeping only 
IPI-diagrams. This is certainly a correct formula for the uni- 
form (i.e. zero momentum) effective action. 

Then one needs the small |u| -expansion of derivatives of R, 
i.e. (|^) as well as 

= i?"'(0+)sign(u) + R""iO+)u + ... (5.29) 
R""{u) = 2R"'{0+)S{u) + i?""(0+) + . . . . (5.30) 

Let us start from: 

E Ri^ab + vf)R{u,d + vf)R{u,f + v^f) . 

abcdef 

We expand in v and of course in diagrams A one must handle 
terms involving R"'{0) and in diagrams B terms proportional 
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to R""{Q). Let us start with diagrams A, which come from 
the following term in the Taylor expansion: 

abcdef 

(5.31) 

Here and in the following, we will drop all combinatorial fac- 
tors. Note that the expectation- values vanish at coinciding 
replicas so there is no need to specify the values of R"'{uab) 
at a = 6. Let us perform the first xy contraction 

R!"MR!"{uaa)R!'M{{vtf{vff{vlff) . 

abcef 

If we now perform a second xy contraction there is a 6aa term 
which is a sloop and thus should be discarded. The 6ad + 
6ba terms build saturated vertices. However the corresponding 
expectation values contain 



i?"'Kd)(«)...)U^a=0 



(5.32) 



which is reasonably set to zero. Thus the first two contractions 
have been performed with no ambiguity leading to 

J2 R"\uab)R!"{Uab)R"M {{vf){vf)iv''/f) (5.33) 
abef 

This term is no more ambiguous. Expanding as in ( I5.29> the 
potentially ambiguous part is 



R"'iO- 



abef 



R!'M{{vf){vf){vlff) 



(5.34) 



clearly free of any ambiguity. It yields the result ( I5.23> . The 
question arises, whether the result may depend on the order. 
We found that when first contracting xy and xz, one repro- 
duces the result ( I5.23t . However when one first contracts xy 
and yz (in any order) one encounters a problem, if one wants 
to contract yz again. The intermediate result after the first two 
contractions is 



E 

abef 



R"'iUab)R"'{Uad)R"{Uaf) (5-35) 



The next contraction between xy contains one term with a sin- 
gle . One would like to argue that this term can be set to 
0. Following this procedure however leads to problems. We 
therefore adopt the rule, that whenever one arrives at a single 
i?"^, one has to stop, and search for a different path. Note 
that this equivalently applies to the recursive constructions 
method. In 2-loop order, one can always find a path, which 
is unambiguous. It seems to fail at 3-loop order; at least we 
have not yet been able to calculate 



(5.36) 




using any other than the sloop elimination method. Whether 
some refinement of the background method can be con- 
structed there is an open question. 



For diagrams of class B one expands as 

J2 R''{uab)R''''{u,<i)R''M{K'f{<)\vl^?) ■ 

abcdef 

Again no need to attribute a value to R""{ucd) for c = d since 
the summand vanishes there. Contract xy: 

J2 R'\u,b)R'''\u,d)R''{u,f){K'){vff{vlf?) ■ 

abdef 

(5.37) 

Contracting yz one gets 

J2R"MR""{Uad) 
abf 

{{vf){vfnR"iUaf)vf R"{udf)vf)) . 

Contracting next xy the danger is the term Sad yielding a sat- 
urated vertex in the middle. But, again if one takes 

R!"'{Uad){{vf)..)\d-.a^Q. (5.38) 

then one gets unambiguously 

Y,R"{Uab)R""{Uab) {{vf){R"{UafW-R"{ubfyj)) . 
ab 

The rest is straightforward. The backgound method thus 
seems to work properly at two loop order 

D. Renormalizability, diagram by diagram 



In section IV A! we have stated that renormahzation diagram 
by diagram gives a method to lift the ambiguity of a given 
diagram, as long as it has sufficient sub-divergences. This 
method is inspired by formal proofs of perturbative renormal- 
izability; the reader may consult 1 9ll S S [H IH 
98] for more details. The key-ingredient is the subtraction op- 
erator S, which acts on the effective action, i.e. all terms gen- 
erated in perturbation-theory, which contribute to the renor- 
malized R, and which subtracts the divergences at a scale /i. 
At 1-loop order, the renormalized disorder Rm at scale m is 
symbolically (with i?o the bare disorder) 



Rn 



where of course the integral <35» depends on m. The oper- 
ator S rewrites this as a function of the renormalized disorder 
R^ at scale /i 



Rni — S 



i?0 + (i?;')'O 



R, + O - O 



Here, the boxed diagram is defined as 



o - o 



(5.39) 



(5.40) 
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The idea behind this construction is that at any order in pertur- 
bation theory, any observable in the renormalized theory can 
be written as perturbative expansion in the bare diagrams, to 
which one applies S. S reorganizes the perturbative expan- 
sion in terms of the renormalized diagrams. The action of 
S is to subtract divergencies, which graphically is denoted by 
drawing a box around each divergent diagram or sub-diagram, 
and to repeat this procedure recursively inside each box. The 
second line of ( I5.39> is manifestly finite, since it contains the 
diagram at scale m minus the diagram at scale /i. This is eas- 



ily interpreted as the 1-loop contribution to the /3-function. 

The power of this method is not revealed before 2-loop or- 
der. Let us give the contribution from the hat-diagram (class 
A): 



Using S, this is rewritten as 




R'oiRo 



Jll\2 



(5.41) 



5R 



(2) 



lll\'2 



i?0 (i?o ) 




R"{R 



lll\2 







(5.42) 



r 



Note that not only the global divergence is subtracted, but also 
the sub-divergence in the bottom loop; and finally the diver- 
gence which remains, after having subtracted the latter (last 
term). Note the factor of 1 = (—1)^ in front of the last dia- 
gram, which comes with the two (nested) boxes. 

Let us halt the discussion of the formal subtraction-operator 
at this point, and not prove that the procedure renders all ex- 
pectation values finite; this task is beyond the scope of this ar- 
ticle, event hough it is not difficult to prove e.g. along the lines 
of once the question of the ambiguities of a diagram are 
settled. However let us discuss, what the subtraction proce- 
dure can contribute to the clarification of the amgibuities. 

In standard field-theory, the main problem to handle is the 
cancellation of divergences, whereas the combinatorics of the 
vertices is usually straightforward. This means that the sum of 
the integrals, represented by the diagrams in the brackets on 
the r.h.s. of ( I5.42> are finite. This ensures of course renormal- 
izability, subject to the condition that all diagrams have the 
same functional dependence on R. Here, the factor R"{R"')^ 
should more completely read 



iR"{u) - R'{0))R"'{u)^ - i^"(u)i?"'(0^ 



(5.43) 



For the first term, there was no problem. However, we have 
seen that the last term was more difficult to obtain. If we 
demand renormalizability diagram by diagram, all diagrams 
have to give the same factor ( I5.43> . Thus, if at least one of 
them can be calculated without ambiguity, we have an unam- 
biguous procedure to calculate all of them. We now demon- 



strate that 



is unambiguous. To this aim, we detail on 



the subtraction operator S, whose action is represented by the 
box. This box tells us to calculate the divergent part of the 
sub-diagram in the box, and to replace everything in the box 
by the counter-term, which here is 



o 



I'i{R"{uf 



2R"{u)R"{Q)) 



(5.44) 



In a second step, one has to calculate the remaining diagram, 
which is obtained by treating the box as a point, i.e. as a local 
vertex. The idea is of course, that the sub-divergence comes 
from parts of the integral, where the distances in the box are 
much smaller than all remaining distances, such that this re- 
placement is justified. Graphically this can be written as 




(5.45) 



Remains to calculate the rightmost term, i.e. to calculate the 
1-loop diagram, from one vertex R{u) and a second vertex 
V{u) R"{uf - 2R"{u)R"{Q). The resuh is in straight- 
forward generalization of 



= h R"{u)V"{u) - R!'{0)V"(u) - i?"(0)T/"(0) 



We need 



V"{u) = R!"{uf + 



(5.46) 



(5.47) 



The omitted terms are proportional to R""R", and contribute 
to class B. We could have avoided their appearance altogether, 
but this would have rendered the notation unnecessarily heavy. 
The term which contributes to ( I5.46t is V"{u) = R"'{u)^. 
It has the same analyticity properties as R"{u), especially 
can unambiguously be continued to u = 0, i.e. V"(Q) = 
R"'{0~^)^. Expression ( I5.46> becomes 



= Ii 



R"{u)R!"{uf-R"(Q)R"'{uf^R"{u)R"'{Q+f 

(5.48) 

without any ambiguity. fl06'] 

To summarize: Using ideas of perturbative renormalizabil- 
ity diagram by diagram, we have been able to compute unam- 
biguously one of the terms in (I5.42> and can use this informa- 
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tion to make the functional dependence of the whole expres- 
sion unambiguous. If we were to chose any other prescription, 
a proof of perturbative renormalizability is doomed to fail, a 
scenario which we vehemently reject. 

E. Recursive construction 

This method is very similar in spirit to the one of section lVDl 
There we had first calculated a subdiagram, and then treated 
the result as a new effective vertex. This procedure can be 
made a prescription, which insures renormalizability and po- 
tentiality, since the 1 -loop diagram insures the latter Only at 
3-loop order appears a new diagram, (I5.36> . which can not be 
handled that way, but the procedure, which is otherwise very 
economic, can handle again most diagrams at 4-loop order, 
using the new 3-loop diagram (I5.36> . 

VI. CORRELATION FUNCTIONS 

Here we address the issue of the calculation of correlation 
functions. We note that it has not been examined in detail in 
previous works on the T = FRG. Usually correlations are 
obtained from tree diagrams using the proper or renormahzed 
vertices from the polynomial expansion of the effective action. 
Thus in a standard theory one could check at this stage that 
correlation functions are rendered finite by the above counter- 
terms, compute them and obtain a universal answer In a more 
conventional theory that would be more or less automatic 

Here, as we point out, it is not so easy. Indeed, as we show 
below if one tries to compute even the simplest 2-point cor- 
relation at non-zero momentum one finds ambiguities already 
at one loop. This is because the effective action (the counter- 
term) is non-analytic. 

Again the requirement of renormalizability and indepen- 
dence of short-scale details guide us toward a proper definition 
of the correlation functions that we can compute. Interest- 
ingly, this definition is very similar as the one obtained from 
an exact solution in the large N limit in f??]. Let us illus- 
trate this on the 2-point function, and at the same time, derive 
the (finite size) scaling function for any elasticity (not done in 
ll67ll ) for massive and finite size scheme. 

A. 2-point function 

We want to compute the 2-point correlation function at T = 0. 
In Fourier-representation it is given by ( 14. U with: 



C{q) = {T('Hq))-,' 



(6.1) 



in terms of the quadratic part of the effective action, which 
reads at any T 

ri^)(^q)^,^l±^Sa, + r^^l{q) (6.2) 



T 

R"{0) 



(6.3) 



replica matrix gives the relation, exact to all orders: 

= -^{R")*{0)m-'^+^'^ . (6.4) 
eh 

R{u) is exactly the function entering the /3-function (in the 
rescaled form R{u)). In the second line we have inserted the 
fixed point form, which thus gives exactly the q — corre- 
lations in the small m limit (i.e. up to subdominant terms in 
1/m) which are bounded because of the small confining mass. 

1 . Calculation of scaling function 

We now compute C{q) for arbitrary but small wave vector 
q, and to one loop, i.e. to next order in e. One expects the 
scaling form ( 14. 2> and that the scaling function is independent 
of the short-scale UV details (i.e. universal), if the theory is 
renormalizable. It satisfies F{0) — 1 and, from scaling should 
satisfy F{z) ^ B / z'^^'^^ at large z. In d = 4 one has F/^{z) — 
1/(1 + z^)^ and we want to obtain the scaling function to the 
next order, i.e. identify 6 in i? = 1 + 6e + O(e^). 

Let us use straight perturbation theory with i?o, defined 
as in Section (IIII C> . including the 1-loop diagrams. This 
amounts to attach two external legs to the 1-loop diagrams 
in Fig.|5] and use a non-anal vtic lllOTll Rq. Our result is: 



{q'+m'YC{q) 



-T^T^^\q 

-R'm - R 



0) 

o'(0 
1 



+ ^2 



(6.5) 



p [p^ + ■m?){{p + qY + w?) 



i.e. by construction here i?"(0) gives the exact off-diagonal 
element of quadratic part of the effective action. Inverting the 



There is however an ambiguity in this calculation, i.e. again 
it is not obvious, a priori, how to interpret the i?o"(0)^ which 
appear. If one computes the one loop correction using \5.2l . 
one must evaluate: 

R"'[Q+f{ulul [ Y.lul-uiWut-ufDGiz-tY . (6.6) 

cd 

One notes that at the very special point z — t there is no 
ambiguity, as the interaction term is analytic to this order 
Then performing the average amounts to take two derivatives 

dn,,du„ of 



which, to this order, is analytic. In this case this is exactly 
the same calculation as for the repeated 1-loop counter-term. 
However, the full expression ( 16. 6> integrated over z^t is, it- 
self, ambiguous. Interestingly, this simple ambiguity already 
to one loop has never been discussed previously. 

Let us first show that renormalizability fixes the form to 
be the one written in ( 16. 5> . Indeed, let us re-express ( 16. 5t 
in terms of the renormalized dimensionless disorder in ( I3.33t 
and JT^ 



<(0) = m^i?"(0) - R"'{{)+Yfr?'I{<d) 



(6.7) 
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As discussed in Section JIII C> . no ambiguity arises when 
taking two derivatives of j3.33> at u — 0+, i.e. the 1-loop 
counter-term is unambiguous. This gives 

= -m' [i?"(0) + R!"[0+fm'{I{q) - /(O)) 

Thus the substitution acts as a counter-term which ex- 
actly subtracts the divergence as it should. The result is finite, 
as required by renormalizability, only with the above choice 
( 16. 5> . Stated differently, the q — {) calculation of ( 16. 5> fixes 
the ambiguity. We show below that the methods described in 
the previous Section also allow to obtain this result unambigu- 
ously. Before that, let us pursue the calculation of the scaling 



function. 

Upon using ( I3.42t and the fixed point equation, we obtain: 



(g2 



1 



1 



(,_2C)^m^(/(g)-/(0)) 
(6.8) 

Apart from the dependence on C,, the calculation of the scaling 
function is very similar to the one given in |67]. We perform 
here a more general calculation which also contains the case 
of elasticity of arbitrary range 



q^ + ^ [q^ + ni^)^ , (6.9) 
and expand in e = 2a — d. Using that, in that case 



-s{p+q/2f-t(p-q/2f-{s+t)m^ 



1 



-^(f)^ Jp Js,t>0 
1 

^ \ 2 J -Jp 



Jp J^yQ 



(1+0 



t q^ 



e/2 



(6.10) 



Defining the 1-loop value of - = (^i + 0(e), we obtain, to the same accuracy, the scaling function in the form (z = \q\/m): 



1 



1 - (1 - 2Ci 



Via) 



1 



e/2 



1 



{l+t)\ 

s)]t-iln(l + s(l-s)z2) 

(6.11) 

(We have used the variable transformation s = 1/(1 + t.)). To obtain b, we need the large z behavior of the scaling function: 



(6.12) 



We want to match at large z: 



y2a 



{l + he + 0{e'))z 



- [l + (e-2C)lnz + 6e + 0(e2)] . (6.13) 



The above result yields 



6 = 6„ = (l-2Ci)[^(|)-V(a)" 

_ / -2(1 - 2Ci)ln2 for a = 1 



-(l-2Ci) 



for a = 2 



(6.14) 



2. Lifting the ambiguity 

Let us now present two additional methods to lift the ambigu- 
ity in the 1 -loop correction to the 2-point function and recover 

In the background method of Section lVCl one performs this 
calculation in presence of a background field, i.e. considering 
that the field u"^ has a uniform background expectation value: 



a a I a 



(6.15) 



with Ua ^ Ub for all a b, and contracting the v^. Then 
at r = the sign of any Ua — ui, is determined, and the 
above ambiguity in ( I6.6I 1 is lifted (contracting further the 
yields extra factors of T and thus is not needed here). Using 
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the background method is physically natural as it amounts to 
compute correlations by adding a small external field which 
splits the degeneracies between ground states whenever they 
occur, as was also found in |72|. On the other hand, perform- 
ing the calculation in the absence of a background field, in 
perturbation theory directly of the non-analytic action yields a 
different result, detailed in Appendix IbI which appears to be 
inconsistent. It presumably only captures correlations within 
a single well. 

The second method is sloop elimination. We want to com- 
pute contractions of: 



(6.16) 



cdef 



where the two disorder are at points z and w respectively. Let 
us restrict to the part proportional to GxzGlw^wy which gives 
the q dependent part of the two point function. Since a is fixed, 
we need to extract the "0-replica part" of the expression af- 
ter contractions (which will necessarily involve excluded ver- 
tices). Starting by contracting twice the two i?'s we get 



c=jtd 



+ ^ R"{ul - ui)R"{ul - ul) 



(6.17) 



Subtracting the sloop W from ( I5.18> gives (up to terms which 
do not depend on both w and z, and which thus disappear after 
the two remaining contractions) 



)R"K 



(6.18) 



Contracting the external u's with ( I6.18> we obtain (restoring 
the correlation-functions) 

/ GuixGwyG'^y 'y^X^abf (6.19) 

The excluded sum can be rewritten as the sum minus the term 
with coinciding indices. Only the latter is a 0-replica term, 
which gives the desired result: 



^i?"'(0+)2 f Gx. 



(6.20) 



This result can also be obtained writing directly the graphs 
with excluded vertices and eliminating the descendants of the 
sloop. 

3. Massless finite size system with periodic boundary 
conditions 

The FRG method described here can also be applied to a 
system of finite size, with e.g. periodic boundary conditions 
m(0) = u{L), and zero mass, which are of interest for numer- 
ical simulations. The momentum integrals in all diagrams are 



then replaced by discrete sums with q = 2im/L, n dlA. One 
must however be careful in specifying the mode q = 0, i.e. 
(tt) — -p- J^Ux- The simplest choice is to constrain (u) = 
in each disorder configuration, which we do for now. Since 
the zero mode is forbidden to fluctuate sums over momentum 
in each internal line exclude q = 0. 

One then finds that the 2-loop FRG equation remains iden- 
tical to ( I3.43> . the only changes being that 

1 . —mdmR has to be replaced by LOlR- 

2. m ^ l/Lin the definition of the rescaled disorder 

3. The 1-loop integral h = j^, ^k^+ra'^yi entering into the 
definition of the rescaled disorder has to be replaced by 
its homologue for periodic boundary conditions: 

1 



used below. 



E 



(27r/L)4(n2 



(6.21) 



Here and below we use a prime to distinguish the different IR 
schemes. As we have seen X is, to dominant order, indepen- 
dent of the IR cutoff procedure. 

Thus we can now compute the 2-point function. Following 
the same procedure as above, we find: 



1 



C{q) = -jL"^ 



l-(e-2C)^(/'(g)-/'(0)) 
with /'(O) i[, and, for q ^ 27m/L: 



E 



1 



(6.22) 



(6.23) 



Thus one finds the finite size scahng function (defined in i4.9\ ) 

c'{d)g,{qL) = q'+'^Giq) = 



2C-e_ 



1 



l-(e-2C)^(/'(g)-/'(0)) 



el- 



(6.24) 



as function of qL = 27rn. The asymptotic behaviour is 
1 



1 



20 ^(/' 



[1 + b'e + (e - 2C) HqL)] 



(6.25) 



which defines b' . The corresponding equation ( I6.13t . when 
regularizing with a mass, holds. Taking the difference be- 
tween the two equations yields 



{b-b')e 



lim 



e-2C 



J'(0) /(O) 



(6.26) 
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where I{q) := I{q)\m=i- To leading order in 1/e, I'{q) = 
/'(O) = I{q) = I{0), such that this difference takes the sim- 
pler form 



2C/ 



(b - b')e = ^^^^ i lim 
e/(0) 



i'{q)-i{q) + m-i'{0) 



(6.27) 



Now observe that for large q the first integral can be bounded 
by 



< 



const 



(6.28) 



which is obtained by estimating the maximal difference of 
integral and discrete sum in each cell (defined by the dis- 
creteness of the sum), and then integrating. The difference 
/(O) — /'(O) is finite and can be evaluated in d = 4 dimen- 
sions. We need the following formulae 



(n^ + 



dne 



nez 



e3,o(e-^) , 



(6.29) 

(6.30) 
(6.31) 



where 63,0 (0 is the elUptic 9-function. This allows to write 
sum and integral as 



E 



1 



(n2)2 



ds s 



(e3,o(e-«))'-i 



(6.32) 
(6.33) 



The difference in question is integrated numerically: 

1 



7(0)-J'(0) 



(e3.o(e-^)) +1 



14.5019 



= -0.00930479 . 



(27r)4 

We thus arrive at (e/(0) = l/iSn^) + 0(e)) 



(6.34) 



6'- 6=™!^ =0.734676(1 -20 (6.35) 

Since the FUG equation and fixed point value {R")* (0) is uni- 
versal to two loops, the final result for the amplitude ratio be- 
tween periodic and massive boundary conditions is 



c{d) e/(0) ^ ' 

= 1 - 1.46935C + 0(e^) 



(6.36) 



B. 4-point functions and higher 

Let us now show how to compute higher correlation-functions 
with no ambiguity using the sloop-method. Let us illustrate 
the method on e.g. the 4-point function: 



(u"(w)t("(a;)w"(y)w"(^))' . 
The following class of diagrams contributes: 



(6.37) 



x,a <I! 



-3^ y,a 



(6.38) 



An arrow indicates contraction towards an external field, with 
position and replica-index as indicated. The combinatorial 
factors are: from the 4 i?'s. the prefactor of the i?'s. 
4! the possibilities to connect the u's to the _R's. 3 for the 
ways to make the loop of i?'s. When contracting first the u's, 
there is another 2^ for the possibilities, to attach the m's to 
the two replicas of R. Therefore only the factor of 3 remains, 
which is the combinatorial factor for ordering 4 points on an 
unoriented ring. 

We start by contracting the four u's with four R's, schemat- 
ically: 



-S>- y,a 



(6.39) 



and then we perform the four other contractions. Since exclu- 
sions at each vertex can be introduced early on, the number of 
possibiUties is not too high and one easily obtains: 

771 r p'« 4 I /I D'« 3 p/// „ p/// 2 p/// 2 . p/// p/// p/// 2 



-L D'" D'" D'" D'" 



(6.40) 



where all terms have to be summed over with excluded repli- 
cas at each vertices. Due to the factors of with an odd 
power, it is not trivial how to project this expression onto the 
space of 0-repUca terms to yield the desired expectation value 
(as in the previous Section a is fixed and thus no free replica 
sum should remain in the final result). 

To perform this projection we will first simplify the above 
expression using sloops. There are a number of possible 
sloops which can be subtracted. The first one is obtained by 
starting from 



(6.41) 
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FIG. 15: Other possible contributions to the connected 4-point func- 
tion. They all are higher replica-terms, and thus do not contribute. 
Arrows indicate contraction towards external points. 



It reads 

The next sloop is 



(6.42) 



(6.43) 



generating 

So 



R. 



ab 



^^ab Rac 



n'" o'" 7?"' ^ 

+ '^^ab^ac'tlad 



The last needed sloop is 

— 



^ab ^ac 
nlll nlll Tjiii nlll — n 



(6.44) 



(6.45) 



leading to 

C d'" 2 nil 

"53 •— -Kq6 



O ??"' R'" ^ R'" _i_ P'" R'" R'" R'" — n 



The simplest combination out of F, Si, S2 and 53 is 



2S2 



3r: 



ab 



(6.46) 



(6.47) 



This expression is unambiguous because only squares of 
R'" (u) appear and it is easily projected onto the 0-replica part 



nlln+\4 



-2R"'{Q 



(6.48) 



e.g. one can replace ^ (1 — Sab)R'ab in (|^^} and use 
free summations. 



Other possible contributions are given on figure [Tsl How- 
ever none of these diagrams contributes. The reason is that 
they are all descendents of a sloop. We start by noting that 



(6.49) 



is a true 1-replica term, i.e. a sloop. When constructing 
a diagram on figure each of the terms in the excluded 
replica formalism is proportional to J6.49I I. thus descendant of 
a sloop. This means, that to any order in perturbation theory, 
at T 0, no diagram contributes to a connected expectation 
value (of a single replica), which has two lines parting from 
one R towards external points. 

Thus the leading contribution in R to the connected 4-point 
function, as determined by the sloop method is the 1-loop di- 
agram 

= — 2i?"'(0^)^ / GrsGstGtuGur 
J rstu 

^ {Gn^rG xsGytG ZU ~^ GujsGxrGytGzU ~1~ Gyjj'GxtGysG zu) 

(6.50) 

If one expressed this result in terms of the force correlator 
i?"'(0+)^ = A"(0+)'* we thus find that this expression is for- 
mally identical to the one that we obtained for the same four 
point function at the T = quasistatic depinning threshold 
(Equation 5.4 in fs?"]). This is quite remarkable given that the 
method of calculation there, i.e. via the non-analytic dynami- 
cal field theory, is very different. Of course the two physical 
situations are different and here one must insert the fixed point 
value for i?"'(0+) from the statics FRG fixed point, while in 
the depinning calculation — A"(0"'') takes a different value at 
the fixed point. In both problems the connected four point 
function starts at order 0(e''). However in some cases the 
difference appears only to the next order in e. For instance, 
we can conclude that the results of L87.] still hold here for the 
static random field to the lowest order in e at which they were 
computed there (of course one expects differences at next or- 
der in e). On the other hand, for the static random bond case, 
the result for the connected four point function will be differ- 
ent from depinning even at leading order in e. It can easily be 
obtained from the above formula following the lines of 187J . 

VII. CONCLUSION 

In this article we constructed the field theory for the statics 
of disordered elastic systems. It is based on the functional 
renormalization group, for which we checked explicitly renor- 
malizability up to two loops. This continuum field theory is 
novel and highly unconventional: Not only is the coupling 
constant a function, but more importantly this function, and 
the resulting effective action of the theory, are non-analytic at 
zero temperature, which requires a non-trivial extension of the 
usual diagrammatic formulation. 

In a first stage, we showed that 2-loop diagrams, and in 
some cases even 1-loop diagrams, are at first sight ambiguous 
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at T = 0. Left unaddressed, this finding by itself puts into 
question all previous studies of the problem. Indeed, nowhere 
in the literature the problem was adressed that even the 1-loop 
corrections to the most basic object in the theory, the 2-point 
function, are naively ambiguous in the T = theory. Since 
the problem is controlled by a zero-temperature fixed point 
there is no way to avoid this issue. An often invoked criticism 
states that the problems are due to the limit of n replicas. 
We would like to point out that even though we use replicas, 
we use them only as a tool in perturbative calculations, which 
could equally well be performed using supersymmetry, or at a 
much heavier cost, using disorder averaged graphs. So repli- 
cas are certainly not at the root of any of the difficulties. In- 
stead the latter originate from the physics of the system, i.e. 
the occasional occurence of quasi-degenerate minima, result- 
ing in ambiguities sensible to the preparation of the system. 
How to deal with this problem within a continuum field the- 
ory is an outstanding issue, and any progress in that direction 
is likely to shed light on other sectors of the theory of disor- 
dered systems and glasses. 

The method we have propose to lift the apparent ambigui- 
ties is based on two constraints: (a) that the theory be renor- 
malizable, i.e. yield universal predictions, and (b) that it re- 
spects the potentiality of the problem, i.e. the fact that all 
forces are derivatives of a potential. Each of these physical 
requirements is sufficient to obtain the /3-function at 2-loop 
order, and the 2-point function and roughness exponent to sec- 
ond order in e. Next, we have proposed several more general, 
more powerful and mutually consistent methods to deal with 
these ambiguous graphs, which work even to higher number 
of loops and allow to compute correlation functions with more 
than two points. We were then able to calculate from our the- 
ory the roughness exponents, as well as some universal am- 
plitudes, for several universality classes to order O(e^). In 
all cases, the predictions improve the agreement with existing 
numerical and exact results, as compared to previous 1-loop 
treatments. We also clarified the situation concerning the uni- 
versality (precise dependence on boundary conditions, inde- 
pendence on small scale details) of various quantities. An- 
other remarkable finding is that the 1-loop contribution to the 
4-point function is formally identical to the one obtained via 
the dynamical calculation at depinning. This hints to a general 
property that all 1-loop diagrams are undistinguishable in the 
statics and at depinning. It would be extremely interesting to 
perform higher precision numerical simulations of the statics, 
and to determine not only exponents but universal amplitudes 
and scaling functions to test the predictions of our theory. We 
strongly encourage such studies. 

Thus in this paper we have proposed an answer to the highly 
non-trivial issue of constructing a renormalizable field theory 
for disordered elastic systems. Contrarily to the closely re- 
lated field theory of depinning, which we were able to build 
from first principles, we have not yet found a first principle 
derivation of the theory for the statics. However, we have 
found that the theory is so highly constrained, and the results 
so encouraging, that we strongly believe that our construction 
of the field is unique. It is after all, often the case in physics 
that the proper field theory is first identified by recurrence to 



higher physical principles such as renormalizability or sym- 
metries, as is exemplified by the Ginsburg-Landau theory of 
superconductivity, for which only later a microscopic deriva- 
tion was found, or gauge theories in particle physics. 

APPENDIX A: SYMMETRIZATION METHOD 

1 . Continuity of the renormalized disorder and 
summary of the method 

The first observation is that one expects (if decomposition in 
p-replica terms is to mean anything) that one can write the 
(local disorder part of the) effective action as a sum over well 
defined p-replica terms in the form: 

p 

= E^ E F(PHua.,.-u.^,, (A.1) 

p ^ ai,...ap 

where the functions F'^p^ have full permutation symmetry. 
The idea of the symmetrization method is that we also expect, 
even at T = that these functions F*^^') should be continuous 
in their arguments when a number of them coincide. 

This seems to be a rather weak and natural assumption. 
Physically these functions can be interpreted as the p-th con- 
nected cumulants of a renormalized disorder, i.e. a random 
potential Vr(u, x) in each environment. Discontinuity of the 
would mean that the Vii{u, x) would not be a continuous 
function. This is not what one expects. Indeed discontinu- 
ity singularities (the shocks) are expected to occur only in the 
force Fii{u, x) — — V„Vr(u, x) as is clear from the study of 
the Burgers equation (see e.g. |54] for the discussion of sim- 
ple case, in the elastic manifold formulation the shocks corre- 
sponds to rare ground state degeneracies). One thus expects 
Vr{u, x) to be a continuous function of u. 

A further and more stringent assumption, discussed above, 
is the absence of supercusp. A supercusp would mean 
i?'(0^) > 0. Thus we assume that the non-analyticity i n the 
effective action starts as The usual interpretation flO^ 
is that there is a finite density of shocks and just counting how 
many shocks there is in a interval between u and u' yields the 
\u — u'\'^ non-analyticity in R{u). 

Let us summarize the method before detailing actual calcu- 
lations. 

We thus define here the symmetrization method assuming 
no supercusp as a working hypothesis. We then compute cor- 
rections to the (local disorder part of the) effective action up 
to a given order in powers of R, with excluded vertices for 
any vector such that Ua ^ Ub for a ^ b, thus with no ambi- 
guity. This yields, as in Section IV B 1! sums over more than 
two replicas with exclusions. These exclusions are not per- 
mutation symmetric so we first rewrite them in an explicitly 
permutation symmetric way which can be done with no ambi- 
guity (see below). Thus we have a sum of terms of the form 

^ rK,,..UaJ, (A.2) 
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where 2 7^ 2 is a short-hand notation for 7^ aj for all i ^ j, 
i.e. symmetrized exclusions. Each function / is fully per- 
mutation symmetric, as indicated by the s superscript. Next 
the non-trivial part is that we explicitly verify that these sym- 
metrized corrections can indeed be continued to coinciding 
points unambiguously, e.g the limit /''(ui, ui, U3, . . . ,Uap) 
exist and is independent of the direction of approach. This 
in itself shows that the continuity discussed above seems to 
work. The existence of a four replica term obliges us to also 
consider three coinciding points. This is done by considering 
/''^(mi, Ui, W3, M4), i.e. symmetrizing the result of two coin- 
ciding points over ui, ^3,^4 and then taking U3 ui. We 
check explicitly that this again gives a function which can be 
continued unambiguously. Thus at first sight, it would appear 
as the ideal method to extract the functions F'^p^ above to or- 
der i?^. 

2. Calculations 

let us reconsider the diagrams of Fig. We first transform 
them in sum with symmetrized constraints. We illustrate this 
on diagram /? where the sum can be reorganized as: 

/3 - ^ KbK'bK'c 

= X! ^abR'ab' + XI ^abR'abR'ac (^-3) 
ab-2^2 abc-2^2 

with clearly no ambiguity. Performing similar rearrangement 
on all the graphs of class A yields the sum of the graphs: 



6aR - 4 X R'abR'ab^ + ^ X! ^abR'abR'a: 
ab:2^2 abc:2^2 

1 nil rtiii T>iii I o 7j" 7j 

^ "-ab^-ac^bc + ^ ^ab^, 



III 2 

ac 



abc;2^2 abc[2^2 

+ 2 X! ^abR'ac^a'd ■ (^-4) 

abcd;2^2 

Now we use the property that has worked on all the examples 
needed here, namely that a symmetric continuous function on 
{{xi, . . . ,Xp);i ^ 3 ^ Xi ^ Xj] is continuous on BP . 
Writing for any /(xi, ..Xp) symmetric and continuous: 

E Vi^^-^a.a,)!{Xa,.-Xa,) (A.5) 

and expanding yields, for the three and four repHca sums: 

^ ^ fabc — ^ ^ fabc ~ 3 ^ ^ faab 2 ^ ^ faaa 
afec,2^2 abc ab a 

^ ^ fabcd — ^ ^ fabcd 6 ^ ^ faabc ~^ 3 ^ ^ faabb 
abcd,2y^2 abed abc ab 



+8 

^ ^ faaab 6 ^ ^ /a 



(A.6) 



in shorthand notations such that fabcd = f{ua, Ub, Uc, Ud) 
This is just combinatorics. 



For the three replica sums the procedure is straightforward, 
as symmetrization makes manifest the continuity. One easily 
finds (we drop an uninteresting single replica term): 



EjyII jylll Tjlll _ 7J" TjlH TjlH Tjll TjlH 2 

^ab^ab^ac ^ ^ab^ab^ac ~ ^ab^ab 

abc;2^2 abc ab 

Ertll JDlll'^ _ \^ pll D'"^ 

^ab^ac / J ^ab^ac 

abc;2^2 abc 



ab 



E 

abc;2^2 



rtll rtlll Tylll _ Till rtlll rtin rtll ((\\ rtll 

^ab^ac^bc — 2^ ^ab^ac'^bc l"7^a 



2 
b I 



abc 



ab 



where in the first line we have applied iA.6\ to fabc = 
sy^abcKbKbKc and so on (we define sym^^ as the sum 
over all permutations divided by p\). 

For the 4-replica term we find that fabcd = 
symabcdR'ab^'acR'ad ^as the following limits (in a sym- 
bolic form, omitting the free summations) 

f — ^ n" (r\\r>'" jy'" ^ n" d'" r?'" 

Jaabc — grt yJj^ab^ac ^ Q^ab^ab^bc 



nil nlll nlll 

^ac^ac^bc 



1 nil nlll^ I ^ nil nlH^ 



faabb — -^R"{0)R'ab 

faaab = -J^RabR"' i^+f + -^R'LbR'ab : 

where at each step we had to symmetrize before taking coin- 
ciding point limits (checking that this limit was unambiguous 
in each case). 

The final result is found to be: 

SARab = «,-i?"(0))«^)2-^<,(i?"'(0+))2) . (A.7) 

The same procedure applied to the repeated counter-term con- 
firms that it is unambiguous and give by ( I3.38> . Thus because 
of the ominous 5/3 coefficient above, rather than the expected 
1 the theory, using this procedure, is not renormalizable. 

Diagrams of class B and C behave properly. One finds with 
the same method their projections on the 2-replica part: 



^/ — ^ 7?"" n" 7?" 

" — 2 ab^ab^ab 

a' = \{2R""{Q)R!'{Q)Rlb + <fc^"(0)') 
7' = /3' 



(A.8) 

(A.9) 
(A. 10) 

5' ^ ^2{R!"'{0)R"{0)R:, + R!2R'LR!'m (A.ll) 

Note the i?""(0) which here is defined as i?""(0) = 
i?""(0+) = i?""(0-) since R!"'{u) can be continued at zero. 
One has, using the expressions given in Appendixlcl 

a" = R!"'mKbf + 2RZR"{Q)Rlb (A.12) 
P" + 6" = -2(i?""(0)i?"(0)<, + RZKbR"m (A.13) 
7" + A" = i?"6i?"(0)' + 2i?""(0)i?"(0X'fc (A.14) 

+ V' = -{R""mKbf + RZiR^'m') (A.15) 



32 



These graphs (more precisely their contribution to 2-rep 
terms) sum exactly to zero: 



a" + /3" + 7" + 5" + ri" + A" + v" = Q 



(A. 16) 



in agreement with the result of the ambiguous diagrammatics 
in the case of an analytic function. 

To conclude, although promising at first sight this method 
is not satisfactory. The projection defined here seems to fail 
to commute with further contractions. For instance one can 
check that upon building diagrams A by contracting the sub- 
diagram (a) in Fig. 12] onto a third vertex does give different 
answers if one first projects (a) or not. Since (a) is the diver- 
gent subdiagram this spoils renormalizability. Since the initial 
assumptions of the method were rather weak and natural, it 
would be interesting to see whether this problem can be better 
understood in order to repak this method. 

APPENDIX B: DIRECT NON-ANALYTIC 
PERTURBATION THEORY 

In this Section we give some details on the method where one 
performs straight perturbation theory using a non-analytic dis- 
order correlator i?o (u) in the action. Expanding in i?o (u), this 
involves computing Gaussian averages of non-analytic func- 
tions, thus we start by giving a short list of formula useful 
for field theory calculation of this Section. One should keep 
in mind that these formula are equally useful for computing 
averages of non-analytic observable in a Gaussian (or more 
generally, analytic) theory. 

1 . Gaussian averages of non-analytic functions: 
formulae 

We start by deriving some auxiUary functions, then give a list 
of expectation values for non-analytic observables of a general 
Gaussian measure. 
We need 



dq (e 



iqx 



iqx\ -7]q 



27J 



(B.l) 



Integrating once over x starting at yields 



1 ^ (.^1- 



1 Jo q 



2arctan ( - ) (B.2) 



The rh.s. reduces in the limit of 77 — > to 7rsgn(a;), which 
gives a representation of sgn(a;) 

sgn(x) = lim — / — sin(ga;)e~'" 

17^0 TT Jo q 

= lim-/ ^sin(gx)e-"l'?l . (B.3) 

By integrating once more, we obtain 

= lim - / ^ [1 - cos(gx)] e""' . (B.4) 
''^0 Ti" Jo 9 



This formula is easily generalized to higher odd powers of |a;|, 
by integrating more often. The result is 



|2ri-l 



2 1"°° dn 

lim-(-l)«r(2n)y^ -^e-''cos(g.)|^ 



(B.5) 

where cos(q2;)|^^ means that one has to subtract the first n 
Taylor-coefficients of cos{qx), such that cos(ga;)|^^ starts at 
order (qcc)^": 



cos(ga;)|^ := ^ 



-{qx) 



2l£ 



(B.6) 



We now study expectation values. We use the measure 



{xx) (xy) 

{yx) (yy) 




(B.7) 



from which the general case can be obtained by simple rescal- 

ing X ^ xj {xx)^^"^, y ^ y/ (yy)^^^- Let us give an ex- 
plicit example (we drop the convergence-generating factor 
e^'''^ since it will turn out to be superfluous.) 



A more interesting example is 

(sgn(a;)sgn(y)) 

Jo q Jo p ,.^±1 



Jo q Jo p ,.^±1 



°° dq ^e-3(p'+'?') (ef«*-e-P«*) 



1 

= ^ / ds / dq I dpe 



-00 1 J-00 P 

t poo poo 



00 J —00 



- ds- 



= — arcsin(t) 



(B.9) 



Another generally valid strategy is to use a path-integral. We 
note the important formula 



ifi^.y)) 



1 



00 00 

X J dx J dyf{x,y) exp 

— 00 —00 



+ — 2txy 
2(l-t2) ' 

(B.IO) 
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An immediate consequence is 

oo 

{f{x)S{y)) = ^^^^1— ^ J dxf{x) exp 



1 

2^ 



2(1-^2) 

— OO 

i=(/(xVl-t2)) (B.ll) 



The very existence of the path-integral representation (IB.10> 
also proves that Wick's theorem remains valid. Let us give 
an example which can be checked by either using ( IB.10> or 

{x''\y\) = {x^){\y\)+2{xy) (a;sgn(y)) 
= {x'){\y\)+2{xyf{6{y)) 



(B.12) 



We finish our excursion by giving a list of useful formulas, 
which can be obtained by both methods: 





2 

TT 


[v/l- 


{xy\y\) 


= 2y 


TT 


{xy\xy\) 


2 

TT 






2 

TT 


(2 + 



VT— + t arcsin(i) 



S^v^l - f2 + (1 + 2*2) arcsin(*) 



(2 + i^) vT^i^ + 3i arcsin(t) 



(B.13) 
(B.14) 
(B.15) 
(B.16) 



2. Perturbative calculation of the 2-point function 
with a non-analytic action 



Let us consider the expansion of the two point function 



in powers of the disorder 

U.EefRoKn, with - 



, ^ (B.17) 

1109II . where 7^ 
ul — u{. We want to 
evaluate these averages at T = with a non-analytic action 
Ro{u). We restrict ourselves to a 7^ 6 since at T = the 
result should be the same for a = b, and we drop the subscript 
from now on. As mentioned above, the Wick theorem still 
applies, thus we can first contract the external legs. The term 
linear in TZ yields the dimensional reduction result ( I2.14t . 
thus we note {u%Uy) = (u°;.u\j)^^ + (u'^Uy) and we find: 

~\G..Gy,, l^R"{uf)R{u'i:^)y^ (B.18) 



up to 0{Ti?) terms. For peace of mind one can introduce the 
restrictions c ^ a,d ^ b m the first sum and c ^ dm the 
second, but this turns out to be immaterial at the end. We 
need only, in addition to ( 15.21 1: 

R'{u) = R"{0)u + ii?"'(0+)u|u| + ii?""(0+)u3 

R"{u) ^ i?"(0) + i?"'(0+)|u| + ii?""(0+)u2 , (B.19) 

6 

since higher order terms in u yield higher powers of T. Using 
( IB.13> to evaluate Gaussian averages this yields: 

{u:uiy^R"'{o+rGix 

/ f GxzGyw y ] 01 {t) — —GxzGyz 



E 

cd 



where we denote: 

G z'in 



t = 
t' 



2Gzz 

G zw 
^Gzz 

2 



{Sab 
(Sac 



Scd ^ Sad — Sbc) 
Sbd — Sad — Sbc) 



01 (i) = -{3t\/l - t2 + (1 + 2t2) arcsin(i) 

TT 

02 (i') = ^ ((2 + - t'2 + M' arcsin(i')) 



Mt') 

(B.20) 

(B.21) 
(B.22) 
(B.23) 
(B.24) 



Note th at th e cross terms _R"(0)i?""(0+) involve analytic av- 
erages iTToll and yield zero (a remnant of dimensional reduc- 
tion). Also, to this order, no terms with negative powers of 
T survive for n = (see discussion below). Performing the 
combinatorics in the replica sums, we find for n = 0: 



«4)' = i?"'(0+)2GLx 



GxzGyw^l 



Gz 



Gz 



GxzGyz^2 



Gz 



$2(5) =-^02(s) + ^02(0 -202(0) 



Gzz 

(B.25) 
(B.26) 

(B.27) 



It is important for the following to note that cancellations oc- 
cur in the small argument behavior of these functions, namely 

one has $i(s) = -s^/tt + O(s^) and $2(5) = s'*/(47r) + 
O(s^). In d = it simplifies (setting Gxy = l/m^ and restor- 
ing the subscript) to: 

/ a b\ -^0(0) A^o"(0^)^ /n/'D3/^12\ ;-T3 oo\ 

{u u ) = -. A 3 \-(J(Rn/m ) (B.2») 

with A ^ (24 - 27\/3 + 87r)/(37r). As such, this formula 
and ( IB.25> seem fine and it may even be possible to check 
them numerically in d = for large m using a bare disorder 
with the proper non-analytic correlator Ro{u). To obtain the 
asymptotic 771 ^ and large scale behavior in any d, one must 
resum higher orders and use an RG procedure. The question 
is whether the above formula ( IB.25> can be used in an RG 
treatment. 
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3. Discussion 

We found that this procedure does not work and we now ex- 
plain why. Let us rewrite the resuh ( IB.25t . including the di- 
mensional reduction term: 



Cabiq) 
Mq) 



-R'^{0)-R'^'iO+nMq)+A,m 
2 



(B.29) 



with i = 1, 2. One notes that if [x) were a constant equal to 
unity, one would recover the result \6.5l obtained in Section 
\y\\ . However, one easily sees that while iIji{x) k, 0.346 ap- 
proaches a constant as a; ^ a where a ^ 1/ A is an ultraviolet 



cutoff, it decreases as "01(2;) 



J2-d 



at large x, as a result 



of the above mentioned cancellations in the small argument 
behavior of the functions ^i{x). Thus the infrared divergence 
responsible for all interesting anomalous dimensions in the 2- 
point function as the non-trivial value of C is killed, and the 
method fails. Even more, the theory would not even be renor- 
malizable. 

We have performed a similar calculation in the dynami- 
cal field theory formulation of the equilibrium problem in the 
limit T 0, using a non-analytic action. There the method 
fails for very similar reasons. Only at the depinning threshold 
we were able to construct the dynamical theory as explained 
in I166. 671 . One might suspect that one has to start with 
a somehow "normal-ordered" theory where self-contractions, 
i.e. terms proportional to Gxx are removed, since they never 
appear in the T — Q perturbation theory. We have not been 
able to find such a formulation. 

Another problem with direct perturbation theory in a non- 
analytic action is that there is a priori no guarantee that it has a 
well defined T = limit. Let us illustrate this on a simple ex- 
ample in d = 0. The following correlation has been computed 
exactly by a completely different method ll69il for the random 
field model in d = (Brownian motion plus quadratic energy 
landscape, (• • • )o indicates averages over all u): 



a result which is also obtained by extrapolation from d = 4 
using the FRG, as detailed in Section lTV B 21 

On the other hand, the above perturbative method yields, 
expanding in a: 



T 



1 



(l + t2) + 0(cr2) (B.31) 



<^ad)/V2- In the zero-temperature limit 
f O(cr^), which is ill behaved. The absence 
of a well defined Taylor expansion in the zero-temperature 
limit is of course a sign that the correct result (IB.30> is simply 



with t 



{5a 



non-analytic in a. Although this solvable example involves a 
correlator Ro{u) with a supercusp, it is possible that a similar 
problem occurs at higher orders (three or higher) in the expan- 
sion of the 2-point function in the case of the usual cusp non- 
analyticity. There have been conflicting claims in the litera- 
ture about this question ll47ll54ll . i.e. the presence of fractional 
powers at higher orders of the expansion in a non-analytic dis- 
order, and it may be worth reexamining. It is however impor- 
tant to note that, since the e-expansion proposed in the main 
text is not based on such a direct expansion, it does not yield 
fractional powers of e, contrarily to what was conjectured in 

Finally, let us point out some properties of non-analytic ob- 
servables. Let us study e.g. Expansion in powers of 

TL yields a first order term ^ 1/\/T. This is the sign of non- 
analytic behavior and indeed it is easy to find that: 



'DR. 



^<(0+)G(0)^ X 

OTT 



\Jl - i2(2 + f) + '6t arcsin(t) - 2 



where ((u^)'^) 



f -Ro(o)'' 



and t 



G{y) 



(B.32) 
The 



first term is obtained by noting that i?o(0) acts as a Gaus- 
sian random force which can then be separated from the non- 
linear force, and the last term, evaluated using the above for- 
mula, is the only one which survives at T = to linear or- 
der in i?o- The formula is interesting as a starting 
point to compute universal ratio, such as {\ua\f'/ ((""a)^) oi" 
{\w^ — w^l)^ / ((Uq — ""b)^)- Indeed one notes that for d < A 
the integral in the term proportional to R'^'{Q^) is infrared di- 
vergent at large y. This is left for future study. 

APPENDIX C: DIAGRAMS OF CLASS C 

In this Appendix we give the expression of each of the dia- 
grams of class C represented in Fig. MA\ in the excluded (non- 
ambiguous) diagrammatics. One finds, including all combina- 
torial factors: 



(B.30) with: 



6 ^ [3 , A = 7 , 



— T)"" T)" Tf" 

" — ~^ab^ac^bc 

13" + 6" = 2i?"J,'i?"fti?ft^ 

_L \" — T>"" T>" D" 
7 + A - -K^fc Rac^ad 



(C.l) 



(C.2) 

(C.3) 
(C.4) 
(C.5) 



APPENDIX D: SLOOP CALCULATION OF 
DIAGRAMS B AND C 

Let us consider the expression 5b R for the B diagrams in the 
excluded diagrammatics ( 15. 14l i. Let us start again from a sin- 
gle sloop ( I5.18> and ( 15.191 1 and contract this time between y 
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and z twice to produce a diagram of type B. This yields: 

•;:;;;;0 -^i?"(o)E<» 

+ X! ^ab^bc + ^ X! ^ab^bc^l 



III! X}ll 

■be 



a^b.b^c 



a^b,b^c 



1 



Ed" p"" td" 



a^b^a^c^c^d 
a^c^b^c^c^d 



^bc^ac ^cd = U 



(D.l) 



the terms R"{0) arise because the first vertex is not contracted 
in the process so one must separate the (unambiguous) diago- 
nal part to obtain excluded sums. 

If one subtracts this identity from ( I5.14t one finds that there 
remain some improper three replica term (the improper four 
replica term however cancels). This is because in the process 
of our last contractions we have generated new sloops, but, 
since replica were excluded they have to be extracted with 
care. 

Let us rewrite the two possible "double sloop" from unre- 
stricted sums to restricted: 

•C^^^;«^^^» — ^ ^ R'ac^ab ^ad 



abed 



a^b,a^c,a^d 



^ac^ab ^ad + 1^ 2^ ^ i^) ^a. 



ini 

b 



a^b 



y^.R"{^)RabKic 
a^b 



(D.2) 



IcO^^ -lE^'XR'bd 

abed 



b 



a^b 



a^b .a^c ,b^d 



a^b 



In the process we have set to zero the terms 

ii?""(0)E<c<'<i->0 (D.4) 

acd 

ii?""(0)E<X.-0, (D.5) 



since they are proper three and four replica terms. 
Defining now: 



<y<:> + 



(D.6) 

The simplest combination which allows to extract the 2- 
replica part is: 



<X> -2.00 

- T^Z^R-ab K^ab - Rq) ■ 



(D.7) 



We now turn to graphs C. The expression for 5cR is given 
as the sum of all contributions in Appendix Icl Within the 
sloop method it gives immediately zero: ScR = 0. This is 
because one can start by contracting the tadpole. Since this is 
a sloop it can be set to zero: 



^G.. E R"Kb)RK,)R{<f)^0 (D.8) 

abcdef 



Upon further contractions, proceeding as in Section fVBl one 
obtains exactly that the sum of all graphs C with excluded 
vertices is identically zero. Graphs C sum to zero since they 
are all descendants of a sloop. 



APPENDIX E: CALCULATION OF AN 
INTEGRAL 

We will illustrate the universality of 

2e{2lA-l!) 



X = 



(E.l) 



using a broad class of IR cutoff functions, namely a propaga- 
tor: 



1 



dx 



9{x) 



+ xm? 



(E.2) 



Here we denote A{x) = / dxg{x)A{x) and we normalize 
/ dxg{x) = 1 (consistent with fixing the elastic coefficient to 
unity). We will show that X = 1 + 0(e) independent of g{x). 
Fkst we write: 
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using the parameterization a = ai + a2 and Xa = (ai — a2)/2 one obtains: 
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The hat-diagram is 

I A - 
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where we split the divergent integral J in pieces, which are either finite or where the divergence can be calculated analytically: 
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This gives the final result for the hat-diagram 
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where we have used that in the 1-loop integral ( IE.4> 
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with a depending on the regulating function g{x) and in the 2-loop integral ( lE.llI 
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with the same a. 
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FIG. 16: Diagrams to order TR^ with excluded vertices. 
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FIG. 17: Diagrams to order TR^ with non-excluded vertices 
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FIG. 19: Divergent unsplitted diagrams to one and two loops. 



APPENDIX F: SUMMARY OF ALL 
NON-AMBIGUOUS DIAGRAMS, FINITE 
TEMPERATURE 

In this Section we give all 1-loop and 2-loop diagrams includ- 
ing finite T, evaluated with the unambiguous diagrammatics, 
which have not been given in the text. We use the unambigu- 
ous vertex J^a^b ^("^ ~ ^b), denote Rat = -R(wa - Wb), 
= R'{ua - Ub) etc.. 

The list of all UV-divergent diagrams up to two loops is 
given in Fig.^] We write their contribution to the effective 
action as 



\ a^b a^b.a^c j 



The total 2-loop contribution is: 



(F.3) 



ah 

The total 1 -loop contribution is 



(F.l) 
(F.2) 



(F.4) 



('21 /2'] 

where 5\' Ris given in ( I5.13> . 5^' Ris given in ( 15.141 1 and 



FIG. 18: Diagrams to order T^R with excluded vertices. 
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while the finite T diagrams are given by 



The repeated counter-term is 
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For an analytic R one substitutes Rab Rab{^ 
above formula and selects the 2-replica terms: 
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Let us show that if one renounces to the projection onto 2- 
replica terms, one can still obtain some formal renormalizabil- 
ity property, but at the cost of introducing an unmanageable 
series of terms with more than two replicas. 

We show how to subtract divergences by adding counter- 
terms of similar form. Let us discuss only T = 0. To cancel 
the 1-loop divergences we introduce the counter-term: 

€^R^iT.l(^abr+ E l^'xAif' (F13) 
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omitting all (excluded) sums. One checks that 

2{5'^g^R + 5fR)^5^^^^^R + 

using that 21 a = if + 0{^). Thus there is some renormal- 
izability property for R. One can thus define formally a /3- 
function: 



(F15) 



mdrnR = eR+ [j2liKbJ 



+S^I^R 
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a^b 



a^b.a^c 



1t2\ 
2h, 
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R however includes a series of terms each with excluded sums 
over p replica. Thus to be consistent one should in principle 
include them from the start and pursue the method. It is not 
clear that it can be closed in any way. 
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[102] One can easily see that the DR result <2.14t arises if one av- 
erages over multiple solutions Ua with some random weights 
Wa ~ I det(V^ + m? + F!u^{x,Ua.)\ (then using the repre- 
sentation of the delta function Gxp(i u{{—V'^ + m^)u + 
F{x,u))) and averaging over disorder using <1.3> . Summing 
over multiple solutions Ua requires instead to include the cru- 
cial weight eicp{~l3H{ua)) in order to select the true ground 
state. 

[103] We thank E. Brezin for a discussion on this point. 

[104] One can also comment on their result for the extremal excur- 
sions A_ff = UinB.^ — Umin. If u Were a Gaussian variable with 
the same 2-point correlator, the exact result for extrema oflog- 
arithmically correlated Gaussian variables predicted in llOdl 

yields: AH = 61nL - cln(lnL) + a where h = 4:,J^, 

c = 27 y'^, 7 = 3/2 and 2 a fluctuating constant of or- 
der 0(1) (in their notations b is our Bd)- Inserting the value 
obtained numerically in l84ll for h yields b = 0.795 and 
c = 0.20. This is in reasonable agreement with the measured 
values 6 ~ .73 quoted in 1 84]. Since deviations from Gaussian 
are not expected to be large, this agreement could probably be 
improved by using the above form of finite size corrections (as 
was done in 1 100] for a one dimensional version where much 
larger sizes had to be considered) rather than the simpler form 
used in [84]. 
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[105] They also measure the decay of the correlation of exp(2i7ru), 
which, within a Gaussian approximation for the distribution 
of It, yields the decay L~^'', with the Bragg glass exponent 

[106] The same procedure can be applied to the dynamics at the de- 
pinning transition. Care has to be taken there, since it exists an 
additional 1-loop counter-term, which is an asymmetric func- 
tion with a vanishing integral. The repeated counter-term at 
2-loop order (integrated over all positions) therefore also van- 
ishes; however it gives a non-zero contribution both to classes 
A and B (chains and hat-diagrams), of which the sum van- 
ishes. In order to ensure finiteness diagram by diagram, these 
contributions may not be neglected. This is discussed in 16^. 

[107] A subtle point in that construction is that if one defines Ro per- 
turbatively from Rtoa given order, then Rq is not the original 
bare action (which is analytic); thus there is no contradiction 
in Ro being non analytic. In a sense, introducing Ro is just a 
trick commonly used in field theory to express a closed equa- 
tion for the flow of R to the same order. The (perturbative) 



exact RG method introduced in f 101*1 does that automatically 
without the need to introduce Ro . 

[108] One should be careful in these arguments, and consider the 
precise definition of R{u). Indeed one could argue that if there 
are many small shocks they could build a supercusp. For in- 
stance consider the non-trivial d = limit of the random field 
model, when V{u) has at large u the statistics of a Brown- 
ian motion. Then, in some definitions of a coarse grained dis- 
order, e.g. such as used in [69] where this model was solved 
exactly, Vr (u) is a continuous one dimensional Brownian mo- 
tion, thus with a infinite number of small shocks and indeed a 
supercusp. However, in the present paper, Vr{u) is defined 
from the replicated effective action, and not from the action, 
and should possess - in that case - a weaker singularity. 

[109] These averages are connected but this is not needed here. 

[1 10] R"{0) can always be set formally to zero by a trivial additive 
random force contribution. 
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